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ABSTRACT 

'•.I 

The  solutions  of  Schroedinger  equations 
(1.1)  e^d^w/dz^  +  =0,  q  =  q(z),  e  -+  0 

"nave  well-known  WKB-approximations ,  but  the  coefficients  in  these  differ 
on  the  two  sides  of  a  turning  point.  A  new  method  fo r  connecting  them 
across  such  points  is  developed  to  extend  present  theory  to  a  more 
general  class  of  turning  points,  which  includes  logarithmic  branch 
points  of  q (z) ,  among  many  others.  To  this  end,  a  delicate  contrac¬ 
tion  for  an  integral  equation  differing  from  those  of  Langer  and  Olver 
is  used  to  show  that  Bessel  functions  can  still  approximate  the  solu-  , 
tions  at  a  certain,  small  distance  from  the  irregular  point  of  (1:1),  ( <7 
even  though  not  uniformly  near  it.  A  novel  feature  of  the  analysis  is 
that  the  extreme  variation  of  the  exponential  kernel  is  here  controlled 
even  on  non-progressive  paths.  Connection  is  completed  radially  by 
means  of  the  same  integral  equation. 
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SIGNIFICANCE  AND  EXPLANATION 

This  work  concerns  the  modulation  of  waves  or  oscillating 
systems,  which  pervade  all  the  science  and  engineering  disciDlines. 
Modulation  occurs  when  waves  travel  through  an  inhomogeneous 
material  in  which  the  local  propagation  velocity  differs  from  place 
to  place,  but  the  differences  are  small  over  a  distance  of  only  a 
wavelength  --  a  very  common  case  in  the  sciences  and  engineering. 

The  resulting  change  to  the  waves  is  mostly  gradual,  but  occasion¬ 
ally  drastic,  as  at  a  shadow-boundary,  where  oscillation  turns  into 
decay  and  quiescence  over  just  a  few  wavelengths.  When  this  phe¬ 
nomenon  can  be  analyzed  via  an  ordinary  differential  equation,  such 
a  boundary  is  called  a  turning  point. 

At  first,  only  the  simplest  turning  points  representing  the 
most  typical  shadow  boundaries  were  studied.  But  then  some  phe¬ 
nomena,  such  as  wave  reflection  and  scattering  cross-sections,  came 
to  be  traced  to  hidden  turning  points  that  become  visible  only  when 
real  distance  (or  time)  is  embedded  in  its  complex  plane.  When  the 
material  properties  vary  in  a  general  manner,  (which  can  often  be 
observed  only  incompletely)  the  hidden  turning  points  can  have 
arbitrarily  complex  structure.  The  following  work  extends  the 
basic  mathematical  formula  for  connecting  waves  with  shadow  across 
a  turning-point  boundary  to  a  much  larger  class  of  variations  in 
the  material  properties  than  had  been  accessible  up  to  now. 
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CONNECTION  OF  CLOSE  QUARTERS  TO  GENERALIZED  TURNING  POLNTS 

+ 

J.  F.  Painter  and  R.  E.  Meyer 

1.  Introduction. 

The  best-known  and  most  important  approximate  solution  to  the 
equation 

(1.1)  e2d2w/dz2  +  q2(z}w(z)  =  0  ,  e-0. 

is  the  NKBJ  solution,  sometimes  also  called  the  Liouville-Oreen 
solution, 

(1.2)  w(z)  ~  Aq~1/,2(z)exp[ie"1  /  q(z')dz'] 

+  Bq~i//2(z)exp[-i£-1  f  q(z')dz']  . 

This  solution  has  only  limited  value  around  a  singularity  of  q'/q- 
It  is  valid  in  a  closed  region  R  which  has  no  singularities  or 
zeroes  of  q  and  which  satisfies  a  certain  convexity  property. 
Because  of  this  convexity  property,  such  a  region  R  can  contain 
only  one  "side"  of  a  singularity  or  zero  of  q.  Inpractice,  oneneeds 
to  .now  w(z)  on  both  sides  at  once,  and  that  requires  a  much  more 
careful  study  of  w(z)  than  needed  to  obtain  (1.2). 
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This  convexity  property  is  defined  in  terms  of  the  natural  inde¬ 
pendent  variable  5,  which  appears  in  (1.2): 


6  =  /  q(z')dz'  . 
0 


Consider  R  =  R^.,  the  region  of  validity  of  (1.2),  as  a  set  in  the 
^-plane.  If  two  points  of  R^.  have  the  same  imaginary  part,  then 
the  horizontal  line  joining  them  must  itself  lie  in  R  As  R 
does  not  contain  any  singularity  or  zero  of  q,  this  means  that  it 
cannot  have  any  5  both  to  the  rig^it  and  to  the  left  of  a  singu¬ 
larity  or  zero  of  q  (Figure  1.1).  The  need  for  this  convexity 
property  arises  in  the  proof  of  (1.2)  from  a  need  to  have  the  expo¬ 
nentials  monotone  on  curves  connecting  points  of  R^,  so  that  they 
can  be  bounded  in  terms  of  their  magnitude  at  endpoints. 

The  limitation  on  where  (1.2)  can  be  applied  leads  to  the  connec¬ 
tion  problem:  given  the  values  of  the  coefficients  A,B  for  which 
(1.2)  applies  on  the  left  side  of  z  =  0,  a  singularity  or  zero  of 
q,  find  their  values  on  the  right  side. 

The  most  satisfactory  way  to  solve  this  problem  would  be  to  get 
a  closer  approximation  to  w(z)  than  (1.2)  -  one  that  is  good  at 
both  sides  of  z  =  0  simultaneously.  Langer  (1931,  1932,  1935), 
Riekstins  (1958) ,  and  Olver  (1977)  have  done  exactly  that  for  the 
simplest  case,  where  q(z)  is  approximately  a  real  power  of  z  near 
z  =  0.  Then  the  approximation  (1.2)  by  exponentials  may  be  replaced 
by  an  approximation  by  Hankel  functions .  The  exponential  approximation 


(1.2)  holds  in  a  set  bounded  away  from  z  =  0,  while  the  Hankel 
finction  approximation  holds  arbitrarily  near  z  =  0  as  well. 

When  q(z)  is  not  a  power  of  z,  the  connection  problem  grows 
more  complicated.  Painter  (1979)  has  shown  for  a  logarithmic  q(z) 
and  imaginary  e  that  the  Hankel  functions  will  not  serve  any  better 
than  the  exponentials  of  (1.2).  The  central  connection  methods  which 
work  when  q(z)  is  a  power  of  z  are  not  adequate  in  greater  gener¬ 
ality  because  the  local  approximands  are  much  less  tractable. 

The  problem  is  that  the  singularity  of  q(z)  at  z  =  0  may,  in 
the  general  case,  be  vicious  enough  to  upset  computations  made  near 
z  =  0  .  The  most  natural  way  to  deal  with  this  problem  is  to  run 
away  from  it,  to  solve  the  connection  problem  as  a  "lateral"  connec¬ 
tion  problem,  solving  the  equation  at  |  z  |  =  ®  . 

This  sort  of  method,  however,  is  doomed  from  the  start.  Distance 
from  the  singularity  makes  it  hard  to  use  information  about  it,  but  it 
turns  out  that  the  answer  to  the  connection  problem  is  strongly 
depende  nt  upon  the  nature  of  the  singularity  at  z  =  0  . 

The  route  to  the  solution  of  this  dilemma  lies  in  the  middle  of  the 
road,  between  the  central  and  lateral  approaches .  Equation  (1.1) 
will  be  solved,  approximately,  at  a  carefully  controlled  distance  from 
z  =  0  .  Some  of  the  advantages  of  both  the  central  and  the  lateral 
connection  approaches  will  show  up  well  enough  to  give  an  answer 
to  the  connection  problem . 

The  method  used  will  be  to  derive,  from  the  diii-runtial  equation 
(1.1),  an  integral  equation  for  the  coefficients  A(£),  B(£)  for 
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which  the  approximate  solution  (1.2)  is  exact  at  z  .  That  is 


(1.3)  w(z)  =  A(C)q_1/2(z)exp(i4/£)  +  B(^)q“1/2  (z)  exp  (-i£ /e)  . 

The  integral  equation  will  then  be  written  into  forms  for  which  it  car 

be  partially  solved  by  the  contraction  mapping  theorem,  i.e.,  by 
iterations.  Some  of  these  versions  of  the  integral  equation  give  A(£) 

y 

or  B(c)  in  various  regions;  others  show  that  A(  £)  -  A  (r)  is  small 

V  v 

in  certain  areas,  where  A  is  the  result  of  substituting  z  for 

q(z)  and  can  be  expressed  in  terms  of  Hankel  functions .  Proving 
that  the  iterations  which  solve  the  equations  work,  involves  many 
estimates,  which  depend  especially  upon  careful  choices  of  the  curves 
along  which  various  integrals  are  computed . 

We  will  be  studying  the  connection  problem  for  a  class  of  singu¬ 
larity  at  z  =  0  defined  in  terms  of  a  modulation  function 
1  2 

<c(0  =  jq'(z)/q  (z)  by  «(£)  ~  y/£  near  zero,  where  y  is  any 
real  constant  <  In  terms  of  q,  this  includes  the  case  which 
Langer  and  others  have  studied,  where  q(z)  ~  z  as  z  -*■  0.  More 
important,  it  generalizes  that  case  to  include  logarithmic  "turning 
points"  such  as  q(z)  ~  zv(log  z)u(v  >  -1)  and  many  others.  It  will 
be  proven  that  the  same  connection  formula  which  holds  for  q(z)  ~  z 
also  holds  for  this  general  case;  that  is,  when  the  region  of  interest 
connects  regions  to  left  and  right  above  the  turning  point  (Fig.  1.1b). 
the  coefficients  Ar ,  in  (1.2)  to  the  right  of  zero  depend  unon 
those  on  the  left  as  follows : 
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(1.4) 


Ar=  kt  -  2B^i  sin(Y-)  +  o  (1) 

Bf  =  +  °(1)  • 

A  coal  formula,  in  which  A  and  B  exchange  roles,  applies 
when  the  connection  is  made  below  the  turning  point  (Olver  1974) . 

Many  other  formulations  of  the  WKB  connection  problem  arise  in  prac¬ 
tice,  but  all  are  reducible  straightforwardly  to  (1.4)  or  its  dual. 

Section  2  is  devoted  to  the  formulation  of  the  integral  equation 
which  was  found  capable  of  solving  the  connection  problem  with  more 
generality  than  before.  In  some  ways,  it  is  intermediate  between  the 
simple  WKB  integral  equation  with  exponential  kernel  (Olver  1974)  and 
Langer's,  with  Bessel  kernel.  The  properties  of  the  new  kernel, 
which  is  related  to  the  incomplete  Ganna  function,  are  discussed  in 
Section  4.  Ihe  main  Section  3  outlines  the  method  of  connection 
along  a  seni-circle  of  radius  4(e)  tending  to  zero  almost  as  fast 
as  e;  Section  5  supplies  details  of  the  proofs.  In  practice,  con¬ 
nection  between  points  left  and  right,  but  so  close  to  the  turning 
point,  is  inadequate  and  Section  6  therefore  extends  the  results  to 
arbitrarily  large  distances. 

While  our  proof  of  (1.4)  is  valid  also  for  y  =  0,  the  formula 
is  not  very  informative  in  that  case.  The  reason  is  that  no  informa¬ 
tion  has  then  been  specified  about  y (f) .  Modulation  functions 
behaving  like  (f  log  0 near  zero  have  been  treated  by  Painter 
(1979)  along  lines  paralleling  those  here  reported,  but  involving 
much  added  labor  because  the  kernel  is  more  complicated.  More 
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definite  results  than  (1.4)  were  deduced  (Painter  1979)  for  that  class, 
which  includes,  e. g. ,  q(z)  which  are,  or  approximate,  a  power  of  log 
z  near  z  =  0. 
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2.  Integral  Equations  for  Connection. 

The  main  integral  equations  needed  to  solve  the  connection  prob¬ 
lem  come  from  the  original  differential  equation,  (1.1),  which  we 
take  to  be  valid  for  z  in  an  open  set  Rz  of  complex  numbers.  We 
assume  q(z)  to  be  holomorohic  and  nonzero  in  Rz;  if  q  have  a 
singularity  or  zero,  it  will  lie  on  the  boundary  of  R  .  We  take  c 
to  be  positive. 

The  WKB  solution  suggests  that  we  look  at  (1.1)  in  a  different 
form.  This  approximate  solution  is 

(2.1)  w(z)  ~  Aq‘1/2(z)exp[i£;/e]  +  Bq"1/2(z)exp [-U/e] 

(2.2)  w’(z)  ~  (iA/e)q1/2exp [i£/e]  -  (iB/e)q1/2exn [~i£/t\ 

where  A  and  B  are  constants  and 

rZ  «  « 

(2.3)  |  =  |  (z)  =  J  q(z  )dz  . 

zo 

With  no  loss  of  generality,  we  may  take  the  fixed  point  zq  =  0. 

The  approximation  (2.1),  (2.2)  is  valid  in  closed  bounded  sub¬ 
sets  S  of  R  with  sufficiently  smooth  boundaries,  which  satisfy 
z  z 

a  "horizontal  convexity"  property.  If  we  set  Sr  =  £(S_)  then  S_ 

c,  Z  c, 

is  horizontally  convex  iff  for  any  two  points  of  with  equal 
imaginary  parts,  the  line  joining  them  also  lies  in  S^.  Because 
the  variable  £  shows  up  in  a  more  fundamental  way  than  does  z 
here,  £  is  a  more  natural  independent  variable. 


This  horizontal  convexity  requirement,  moreover,  is  the  source 
of  the  connection  problem.  Suppose  that  q(z)  has  a  possible  singu¬ 
larity  or  zero  at  one  point  of  the  boundary  of  R  ,  say  at  z  =  0. 
(Here  and  henceforth  we  as sene  that  q(z)  is  holomorphic  and  nonzero 
at  the  other  boundary  points  of  R  ).  We  must  assume  q(z)  to  be 
integrable  at  z  =  0  so  f  can  be  defined  by  (2.3).  We  may  use 

(2.1)  to  compute  w  in  a  large  neighborhood  of  a  point  to  the  left 

of  zero,  s;  <0,  and  in  a  neighborhood  of  a  point  to  the  right  of 

zero,  f  '>  0;  nevertheless  if  we  were  to  know  the  coefficients 
r 

A, , B;  for  which  (2.1)  holds  near  ,  equation  (2.1)  would  not  help 

us  find  the  correspondin'  coefficients  A  ,B  for  .  This  is 

r  r  r 

because  (2.1)  cannot  be  applied  in  a  set  which  contains  both  r  and 

at  once . 
r 

For  the  connection  problem  we  seek  to  know  and  in  terms 

of  A„  and  By.  We  will  find  Ar  and  Br  by  considering  (2.1), 

(2.2)  as  equations  exactly  true  throughout  R  ,  that  is, 

q1,/2w(z)  =  A(Oexp(i^/e)  +  B(Oexo(-ii/e) 

(2.4)  ev' (z)/(iq^2)  =  A(C)exp(iF/e)  -  B(Oexo(-i£/e) 

for  (e  Rf  =  f(Rz)  • 

Thus,  A^  =  A(E£).  B^  =  B(^),  Ar  =  A(fr) ,  and  Br  =  B^p . 

Substitution  of  (2.4)  into  (1.1)  yields  equations  for  the  modula¬ 
tion  coefficients  A(£)  and  B(f) , 
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(2.5)  dA/d4  =  ?(£)B(i)e"p* 

(2*6)  dB/d4  =  ?(|)A(£)ep* 

where  <p  is  the  "modulation  function" 

(2.7)  <p  =  ■~q_2dq/dz  =  q_1dq/d« 

and  p  =  2i/e  . 

In  order  to  solve  (2.5)  and  (2.6),  we  shall  associate  them 
with  integral  equations.  First,  though,  we  must  ensure  that  the 
integral  equations  will  exist:  we  assume  that  Rr  is  connected 
and  simply  connected,  and  that  it  contains  the  set 

:  |  Re|  |  s  5q,  Im£  s  5Q}  \  {%  :Re£  =  0,  0  >  Im£  >  -  ®} 

with  6C  >  0  (Figure  2.1).  For  example,  if  q(z)  be  singular  only 
at  z  =  0,  could  consist  of  all  complex  numbers,  less  z  =  0 

and  a  branch  cut  from  zero  to  infinity.  We  also  assume  that  ?(£,) 
grows  no  more  than  exponentially  as  |  — •  -  i  »  for  -  6qs  Re|  s  6q 
(Re|  /  0)  •  That  is,  we  assume  that  ?>(|)  =  ®  (exp  0c(|  I  ))  with 

k  a  constant  (not  depending  on  £)• 

At  this  point,  a  choice  has  been  made  in  regard  to  the  cut 


(Fig.  2.1),  and  its  implication  should  be  explained.  First  of  all, 
the  original  formulation  of  the  Schroedinger  equation  (1.1)  is  some¬ 
what  arbitrary.  It  is  the  most  general,  linear  ordinary  differential 


equation  of  second  order  in  normal  form,  but  there  is  a  whole  family 
of  "Liouville"  -  transformations  z  -*■  z'(z)  leaving  the  normal  form 
of  (1.1)  invariant,  but  of  course,  changing  the  coefficient  q(z) 
into  a  different  coefficient  Q(z')  and  changing  the  appearance  of 
a  connection  problem.  This  explains  why  connection  problems  can  come 
in  so  many  guises  and  why  it  is  difficult  to  tell  at  the  outset 
whether  two  people  are  talking  about  (Liouville  transforms  of)  the 
same  connection  problem  or  about  genuinely  different  ones. 

By  contrast  to  z,  £(z)  is  the  distinguished  variable  charac¬ 
terizing  the  oscillatory  nature  of  the  solutions.  Where  an  approxi¬ 
mation  (2.1),  (2.2)  holds,  |A|  and  |B|  represent  amplitudes  and 
-1/2  -1/2) 

-e  arg(Aq  )  and  -e  arg(Bq  ,  phases;  £/2tte)  measures 
distance  or  time  in  local  wavelengths  or  periods.  Any  canonical 
formulation  of  the  connection  problem  must  therefore  by  in  terms  of 
£.  The  real  axis  of  £  is  the  line  on  which  (2.1),  (2.2)  are  purely 
oscillatory,  without  exponential  growth  or  decay,  and  this  makes 
connection  between  <  0  £^  >  0  the  nest  conum  form  of  the 

problem.  It  leaves  a  pair  of  problems,  however,  one  with  the  cut  in 
the  lower  half -plane  and  the  other,  with  the  cut  in  the  upper  one. 
These  are  analogous,  dual  problems,  differing  just  by  an  exchange 
in  the  roles  of  A(£)  and  B(£).  For  definiteness,  that  with  the 
cut  in  the  lower  half-plane  is  selected  for  study  here. 

The  integral  equations 

(2.8)  B(£)  -  Bj  =  ®(s)A(s)epSds 

h 
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(2.9) 


A(|)-A£  =  <p( s)B(s)e'PSds 

^  l 

are  imiEdiately  associated  with  (2.6)  and  (2.5).  For  the  present 
choice  of  cut,  however,  a  single  integral  equation  involving  A 
alone,  rather  than  A  and  B  as  m  (2.9),  is  more  helpful.  Towards 
that  end,  integrate  (2.9)  by  parts: 


s  s  =  |  £  s 

(2.10)  A(|)-A,  =  B(s)  /  Wt)e‘Ptdt|  -  J  ^§lds  /  9»(t)e-ptcit 

4  s  =  S*  4£  I 

To  clarify  the  structure  of  (2.10),  let 

(2.11)  J(4)e_P|  =  f  e”Pto(t)dt 

I 

For  sufficiently  small  e,  the  exponential  factor  decreases  much 
faster  than  *  increases  as  t-*--i»,  so  that  this  integral  (2.11) 
converges.  The  kernel  in  (2.10)  is  then 

/ 

(2.12)  )(|,s)e"Ps  =  /^e"pt^(t)dt=  j(s)e"ps  -  )(|)e‘p^ 

s 

and  by  (2.6),  (2.10)  becomes 

(2.13)  A(|)  -  j(|,4£)exp(-p4^)  +  /  p(s)A(s)j(4  ,  s)ds 
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This  will  be  our  basic  integral  equation.  We  shall  solve  it  for 
A(" ),  returning  to  the  much  simpler  problem  of  finding  B(rr)  only 
at  the  end. 

We  will  usually  take  the  path  of  integration  for  (2.8)  and  (2.13) 
to  be  an  arc  of  the  circle  js|  =  5.  The  radius  6  =  ' (>:)  will  be 
chosen  as  a  compromise  between  two  conflicting  needs.  First,  if  q 
has  a  singularity  or  zero  at  =  0,  that  will  determine  the  connec¬ 
tion  coefficients,  so  we  must  do  computations  near  zero  in  order  to 
get  an  answer.  On  the  other  hand,  if  q  has  a  singularity  or  zero 
at  E,  =  0,  then  y  nay  be  behaved  badly  enough  to  make  computa¬ 
tions  excessively  difficult  if  they  be  done  too  near  r  =  0.  Such 
conflicts  will  lie  at  the  heart  of  much  of  the  analysis. 

The  basic  conditions  on  5(c)  are  that: 

(2.14)  (a)  6(e)  >  0  (b)  5(e)  -  0  (c)  6(e)/e  -  -  as  e-o 

Condition  (2.14b)  is  strengthened  in  Section  3,  equation  3.21. 

We  take  and  %x  to  lie  on  the  circle  |||  =  5;  i-e., 

(2.15)  C£=-5  yr=S 

The  coefficients  A(|),  B(|)  at  %  =  |£  and  %  =  %x  are  normally 
approximately  equal  to  the  coefficients  at  greater  distances  from  zero, 
such  a.e  |  =  -1  and  |  =  1.  That  will  be  discussed  in  Section  6. 
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3. 


Close  Connection. 


3.1.  In  order  to  approximate  the  WKB  coefficients  on  the  right, 
Ar  and  B^,  in  terms  of  given  values  of  the  coefficients  on  the 
left,  and  B^,  we  must  knew  something  about  the  singularity 
between  the  right  and  left,  at  £  =  0.  Since  we  are  interested  in 
q(z)  which  have  logarithmic  behavior  at  z  =  0,  we  choose  an 
assumption  which  will  cover  such  q.  We  assume  that: 

(A)  fyfiiO  -*•  Y  as  1 5 1  -»•  0 

uniformly  in  arg  £  for  which  *  is  defined;  i.e.,  E,  e  R,_  Here 
y  is  an  arbitrary  constant,  with  -°°  <  y  <  j. 

In  Section  5.1,  it  wall  be  shown  that  (A)  holds  wdth 
y  =  j  v(v  +  1)  ^  when  q(z)  =  zv(log  z)u,  wzith  v  and  y  real 
and  v  >  -1.  (A)  also  covers  the  "fractional"  turning  point  class 

studied  by  Langer,  Riekstins  and  Olver  (1977). 

To  use  the  basic  integral  equation  (2.13),  we  must  begin  by 
estimating  the  kernel  j(£,s),  which  is  complicated  by  the  lack  of 
information  on  the  rate  of  approach  to  the  limit  (A) ;  the  full  report 
is  postponed  to  Section  4.  The  method  is  to  substitute  (A)  into 
(2.11)  and  (2.12).  Once  the  resulting  error  terms  are  successfully 
dealth  with,  j(?)e-p^  will  be  approximated  by  an  incomplete  gamma 
function: 


-pt 


(Y/t)dt 


This  can  be  further  approximated  by  a  well-known  formula,  except 
near  £  =  .  At  £  =  £^,  another  method  can  be  used.  The  most 

important  difficulty  in  dealing  with  the  error  terms  involved  lies  in 
the  fact  that  exp(-  pt)  in  (2.11)  is  large  when  the  imagainary  part 
of  t  does  not  quickly  approach  zero;  while  if  t  does  quickly 
approach  zero,  then  not  only  will  ,p(t)  in  (2.11)  grow  large,  but 
so  will  the  error  term  <p(t)  -  (Y/t)  =  o(t“*).  This  sort  of  problem  is 
handled  primarily  by  care  in  choosing  the  different  paths  of  integration 
on  whicn  (2.11),  (2.12)  and  their  error  terms  are  computed.  The 
results  are  summarized  below.  These  asymptotic  formulae  hold 
uniformly  as  e  **  0,  for  £  and  s  lying  on  arcs  of  the  circle 
|£  I  =  1  s |  =  6(e);  that  is,  1  PI  I  =  I  p  s|  =  26/e  .  These  arcs  are 
pictured  in  Figure  3.1. 

(3-1)  J(|)=  (Y  +  o(l))/(p|  )  for  Os  arg£  s  tt  -  ©o (Figure  3.1  (a)) 

where  0Q  >  0  is  an  arbitrary  constant. 

(3.2)  J(S,s)e"ps=  (Y  +  o(l))[e'ps/(ps)  -  e"PV(P|)] 

for  0  s  arg£  <  it  -  0Q  and  Osargs<ir-©o  (Figure  3.1  (a))  . 

-P&# 

(3-3)  j(£,£j)e  =  (Y  +  o(l))[exp(-p|r)/(p|r)  -  exp(-p|)/(p|)-2ri] 

for  0«  arg£  S  tt  -  ©Q  (Figure  3.1  (a)) 

(3.4)  j{£,  s)e  PS  =  ®(e”PV(p|)) 

for  j  s  arg£  s  arg  s  s  w  (Figures  3.1  (b)  and  (c))  . 
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3.2.  We  need  a  bound  on  A(r)  before  we  can  compute 
it.  First,  note  that  the  derivation  of  (2.13)  can  be  repeated  with 
£r  in  place  of  ^  to  give 

(3.5)  A($)  -  Af=  Brj(4,Sr)exp(-p$r)  +  <p(s)A(s)j(|  ,  s)ds 

There  is  exponential  behavior  in  A;  we  extract  most  of  it 
by  multiplying  A(0  by  an  exponentially  snail  factor: 

(3.6)  a(£ )  =  A(£)eP^  . 


Then  equations  (2.13)  and  (3.5)  become: 


(3.7)  a(|)  -  adexp(p|  -  Pld)  =  Bd  ]{%  ,  |d)  exp(p|  -  p|d> 


a(s)j(| ,  s)exp(p|  -  p  s)  <p(s)ds 


where  "d"  can  be  either  "  P "  or  "r".  To  have  the 
exponential  in  the  integral  in  (3.7)  bounded,  we  require  £  to 
lie  in  the  same  quadrant  as  |d  (Figure  3.2) .  That  is,  if  |d  = 
then  ?i  arg|  <  it,  so  that  it  s  arg(p|)  i  arg(Ps)  s  arg(p£  ,)  =  3^/2. 

*  h 

If  !d  =  Sr.  ^en  0  s  arg|  «  J,  so  that  arg(p£r)  =  tt/2  *  arg(Ps) 

<  arg(p£)  £  tt  . 

With  such  £,  if  s  lies  between  %  and  *d  as  in  the 
integral  of  (3.7),  then  ps  will  have  greater  real  part  than  p| , 


so  that  exp(p|  -  ps)  will  be  bounded.  Moreover,  we  may  use  (3.2) 


or  (3.4)  to  find  that,  for  s  in  the  integral  in  (3.7)  or  for  s  = 

(3.8)  j(4,  s)exp(p|  -  ps)  =  ©(p^f1  =  ©(e/6)  -  0 
Substitute  (3.8)  and  assumption  (A)  into  (3.7). 

(3.9)  a(e)  =  +  Bd©{£/5)  +  |  ia||f  ®(|  /  ^) 

where  ||a||^  is  the  sud  norm,  applied  for  £  or  s  on  the  quarter- 
circle  allowed  (Figure  3.2);  that  is,  ||a||^  =  sup{|a(s)|  ;  j  <  arg(s) 

<  tt,  |s  |  =  6}  and  ||a||^  =  sup{  |a(s)  |  :  0  £  arg(s)  <  tt/2,  |s  |  =  6}. 

Let 

(3.10)  C  =  max(|A£|,|ArJ,|B£|,|Br|)  . 

We  shall  see,  of  course,  that  C  is  usually  bounded  as  e  -*•  0, 
but  that  is  not  yet  known.  The  integral  in  (3.9)  is  bounded  by  tt/2, 
and  the  exponential  decreases  (or  is  constant)  in  magnitude  as  z  -*■  0, 
so  that  (3.9)  becomes 


a(|)  =  ©(C)  +  ©f|)  |  |a|  |f 


Therefore , 
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I  I a| if  =  ©(C) 

and,  hence, 

(3.11)  Mai  1^=  sup  [\a(i)\  :!§|  =  6,0s  arg|  s  n}  =  ©(C) 

We  may  use  (3.6)  to  translate  this  back  into  the  original  A(£) 
notation . 

(3.12)  A(£)  =  C  ©(e~pS)  =  C  ©(e26/£) 

for  |  on  the  entire  semicircle  ||  |  =  6,  Os  arg£  s  ir 
(Figure  3.2(c))  . 

The  argument  that  produced  this  bound  can  be  carried  out  more 
precisely  to  give  a  formula  for  A(0  valid  on  the  arc 
2  it  _<  arg  Cl  11  ~  9q  (Figure  3.1(d)).  That  is  done  in  Section  5.3. 

3.3  Now  we  can  compute  A(C)  by  showing  it  to  be  close  to  a 
knovn  function  AY(C) .  We  have  described  A(0  by  following  three 
equations : 

,6 

(2.13)  A(| )  =  +  B.  J(|,e^)exp(-P|i)  +  £  A(s)j($,  sMs)ds 
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(2.12) 


j(4,  s)exp(-ps)  =  /  exp(-pt  )<?(t)dt 
s 

(A)  <p(  t)  =  (V  +  o(l))/t 

We  define  AY(£ )  byreplacing  ,?(t)  by  exactly  Y/t ;  that  is, 

(3.13)  AY(*  )  =  A.  +  B.  jY(4,|  ,)exp(-pt,)-L  /  AY(s)j\,s)Y  s  *ds 

Z  L  *•  £ 

v  t 

(3.14)  j  (4  ,  s)exp  (-ps)  =  J  exp(-pt)Yt  dt 

s 

If  q  (Z)=ZV,  then  p(t)  =  Y/t  with  Y  =  |v(l+vf1;  thus 
AY(£)  is  the  modulation  coefficient  for  the  special  case  of 
q(z)  =  zv  .  For  such  q,  the  solution  w(z)  to  (2.1)  can  be 

written  exolicitly  in  terms  of  Hankel  functions  (Section  5.2),  so 
we  know  exactly. 

In  Section  4.8,  we  show  that  j  and  j  '  are  approximately 
equal  in  the  following  sense.  The  error  bomds  are  defined  in 
terms  of 

(3.15)  h(e)  =  max(g(6),  exp(-2S/e  )) 

where  g(6)  is  the  maximum  error  in  assumption  (A),  that  is, 

(3.16)  g(5)  =  sup  f  \Z  <p(t)  -  Y  !  :  !|  I  S  5,  £  e  } 


T" 


r 
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Note  that  g(6)  -  0  as  6-0,  so  also  h(e)  -  0  as  £  -  0  • 
All  of  the  following  3  equations  hold  for  |  and  s  satisfying 
||  (  =  |  s |  =  6,  and  with  e  —  0  . 


(3.17)  j(£,  s)  =  (1  +  &(h(s)))jY(| ,  s) 

for  0  s  arg£  s  t  -  0  ,  0  £  arg  s  £  it  -  e_  (Figure  3.1  (a)) 

3  o  o 

(3.18)  =  (l  +  &(h»iY(!.!£  ) 
for  OS  arg  £  *  ir  -  eQ  (Figure  3.1  (a)) 

(3-19)  j(|,  s)  =  jY(|,s)  +  ®(h)  «  6-1exp  (ps  -  p§  ) 

for  -|t  s  arg|  sc  arg  s  s  tt  (Figure  3.1  (c)). 

y 

Since  we  want  to  use  this  to  show  that  A(£)  -  A  (£)  is  small, 
we  next  combine  (2.13)  and  (3.13)  into  an  equation  for 
A(|)-AV(S). 

(3.20)  A(|)  -  AV(t)  =  B/  ~  JV(€.S|)]e*P(-pSi)  + 

+  f  AV(s)[j(4  ,  s)<p(s)  -  jV(|(  s)ys-1]ds  + 

ii 

+  f1  [A(s)  -  AY(s)]j(|,  s)<p(s)ds 

%  i 

This  equation  will  be  solved  approximately  in  the  next  two  subsections. 
They  will  depend  upon  a  restriction  on  the  choice  of  the  arbitrary 
function  5(e).  So  far  we  have  required  that  £  <  5(e)  <  1  .  Nov. 
we  also  require  5(£)  to  satisfy 


(3-21) 


£5~*g(6)  exp  (26/s  )  —  0  as  £  —  0 

Consequently,  h{£)  also  satisfies  (3.21).  Whatever  the  function 
g(6)  may  be,  this  condition  (3.21)  can  be  satisfied  by  some  choice 
of  6(e).  For  example,  we  could  define  a  function  e(5)  by 

(3.22)  £(6)=  26 1  log(g(5))|  _1 

This  decreases  nrmotonically  to  zero  as  6  0  because  g  ( 5 )  does; 

thus  e(6)  can  be  inverted  to  give  6(£)  satisfying  both 
e  <  6  <  1  and  (3.21). 

In  the  case  of  greatest  interest,  where  q(z)  =  zv  (log  z)^,  we 
have  g  (6 )  =  &(log  6)”1,  by  (5.7),  so  that  the  requirements  for 
6(e)  can  be  satisfied  with 

6e  log  |  log  6 1 

This  increases  only  slowly;  i.e.,  equation  (3.21)  has  forced  us  to 
put  6(t  )  rather  near  e  . 

3. A.  In  this  section  we  shall  solve  (  3.20)  for  %  in  the  arc 
III  =  6,  -|ir  <  arg£  <  it  (Figure  3.3).  To  bound  the  first  bracketed 

term  in  (3.20),  we  use  equation  (3.19)  with  s  =  | 

l 

(3.23)  [j(|  )  -  jY(| ,  |^3exp(-  p% ^  )  =  <Mh)  «  5-1exp  (-p|)  . 
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The  integrand  in  the  second  terra  of  (3.20)  contains  ,  which  is 
bounded  by  (3.12).  For  the  case  q(z)  =  zv  which  leads  to  A’, 
it  is  already  known,  as  will  be  discussed  in  section  5.2,  that 
equation  (2.2)  (with  A  and  B  constants)  holds  as  close  to 
zero  as  |||  =  6,  and  that  the  connection  formula  (1.4)  holds; 
therefore,  C  defined  by  (3.10)  satisfies  C=  ©(1)  as  £  —  0, 
and  (3.12)  reduces  to 

(3.24)  AV(£)  =  ©(exp(-p|)) 

Substitute  (A),  (3-24),  (3.19),  and  (3.2)  into  the  second  term  cf 
(3.20). 

J  AY(s)[j(£,  sMs)  -  jV(£,  s)Y  s-1]  ds  = 

h 

^  V  —1 

=  f  ®(exp(- ps))[©'g(5  ));(£,  s)  +  j(£,s)  -  j  (|,  s)]Ys_  ds 

(3.25)  ^ 

=  +  ©(h)  €  5"1e”P$]Ys”1ds 

h 

=  ©(h)  <  5_1e'p*  . 

The  last  term  of  (3-20)  can  be  estimated  by  the  same  method  a 
in  section  2,  again  using  equation  (3.S)  .  Define 

(3.26)  b(| )  =  [A(| )  -  a\| )]  exp  (p£ ) 
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then  (3.20),  (3.23),  and  (3.25)  may  be  combined  to  give 

-i  £> 

(3.27)  b(£)  =  ©(h)  £  5  +  /  b(s)exp(p|  -  ps)j(|,s)<p(s)ds 

Now  substitute  (3.8),  and  use  the  norm  definition  after  (3.9)  and 
the  fact  that  J  <p( s)ds  =  ©(1)  f  s"*ds  =  ©(1). 

-1  ^  \ 

(3.28)  b<£)  =  ©(h)  c  6  +  |  |b|  \m  ©(E6'1) 

It  follows  that 


(3.29)  b(£ )  =  ©(h)  €  /6 


By  (3.26),  this  may  be  translated  back  into  our  original  terminology . 
For  the  |  allowed  in  this  section  (Figure  3.3), 


Now  equation  (3.21)  gives  us  the  result,  that  A^(£)  is  a  first 
approximation  to  A(|).  For  %  in  the  upper- left  quarter-circle 
(Figure  3.3): 

(3.31)  A(|)=  AY(£)  +  o(l) 
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3.5.  In  this  section  we  shall  solve  (3.20)  for  |  in  the  arc 
III  =6,  0  i  arg |  t  j -n  (Figure  3.3).  That  will  give  us  the  solution 
to  the  connection  problem.  The  first  bracketed  term  in  (3.20)  may  be 
bounded  by  means  of  equations  (3.18),  (3.3),  and  (3.21). 

[j(6,5x)  -  JY(^e^)]  exp  ( -  P  )  =  fe(h)j(|,|i)exp(-p|^) 

=  ®(h)[®(l)  +  S(p|)_1exp(-  p|)] 

(3.32) 

=  ®(h)[®(l)  +  ®(e  6"1  exp  (25/e  ))] 

=  o(l) 

For  the  second  term  of  (3.20)  we  start  by  substituting  assumption 
(A)  and  equations  (3.16)  and  (3.24) 

f*Ay(sm.sM')  -  s)Ys_1]ds 

(3.33)  1 

=  /  ®(e"PS)[$(g)j(4,  s)  +  J(|,  s)-jY(|,  s)]YS"lds 

£ l 

By  (3.2)  and  (3.4),  we  can  take  care  of  the  first  term  in  brackets: 

(3.34)  j(|,  s)e"ps  =  ®{(p|f1e"p^)  +  ® ((Ps)_1e"ps) 

=  ^(eS”1  exp  (26/s)) 


For  the  rest  of  the  bracketed  expression,  we  deal  separately  with  s 
in  each  of  the  two  quadrants  in  Figure  3.3. 


If  s  be  In  the  right  quadrant,  0  *  arg  s  a  ^  »  then  by  (3.17) 

and  (3.2) 

[j(4.  s)  -  jY(|,s)]e"ps  =  ©(h)j(g,  s)e'Ps=  ©(h)(psf 1  e“ps 

(3*35)  _i 

=  ©(h)  «  6  exp  (26/e) 

On  the  other  hand,  if  s  be  in  the  left  quadrant,  -^ir  a  s  a  it  ,  then 
we  must  first  split  up  j(|,  s),  before  applying  equations  (3.18), 
(3.3),  and  (3.19). 

[)(£,  s)  -  jY(|,  s)]e_Ps  =  [!(£,£,£  )exp(-P^)  -  j(s,S£)exp(-p££)j 

-  DY(^l£)exp(-p^£)  -  jY(s,^)exp(-p^;)] 
(3.36)  =  fi(h)J(|  $^)exp(-p£j^  )  -  ©(h)  «  6_1exp(-ps) 

=  ©(h)[e6~1exp(-p^)  +  ©(1)]  -  ©(h)  e  6_1exp  (-  ps) 

=  © (h)[«  5"1exp(26/e)  +  ©(1)] 

Now  we  can  put  together  equations  (3.34),  (3.35)  and  (3.36),  by 
substituting  them  into  (3.33).  Then  use  (3.21). 

f1  AY(s)[j(£,  sMs)  -  jY(|,  s)Vs_1Jds 

*1  . 

=  f  ©(h)  t  6~*  exp  (25/e }s  *ds 

u 

=  0(1) 


Next  break  up  the  last  integral  in  (3.17)  into  two  parts,  so  that  when 
(3.32)  and  (3.36)  are  substituted,  it  becomes: 


V  t*  V 

(3.37)  [A(£)-  AY(4)]  =  o(l)  +  f  [A(s)  -  AY(s)]j(|,  s)^(s)ds 

A  V 

+  J  [A(s)  -  Ay(s)]  j(£  s)(p(s)ds 
i  6 

The  first  integral  in  (3.37)  can  be  handled  with  (3.2)  and  (3.30), 
the  result  of  the  previous  section. 

16  Y  16 

f  [A(s)  -  Ar(s)]j(|,  sMs)ds  =  f  <&(h(e))(£/&)e”  "sj(£,s)<?(s)  c's 
4  ^  4£ 

=  S(h(e))(£/5)  A Y  +  odn^psfVP3  -  (P4 f  )s^ds 

(3.38)  1 

=  S(h(e))(£/5)2exp(26/E) 

=  o(e/5) 

\ 

In  the  last  integral  in  (3.37)  we  may  estimate  j(£,s)  by  (3.2)  and 
the  fact  that  exn(-pO  is  smaller  than  exo(-ps)  there,  and  we  may 
estimate  <p(s)  by  (A).  The  result  is  that  (3.37)  and  (3.38)  imply 

(3.39)  [A(|)  -  AV(4)]  =  o(l)  +  ©(e/6)l|A(s)  -  AY(s)1|w 

It  follows  that  in  the  quadrant  studied  in  this  section,  and  hence  in 


the  entire  semicircle  |||  =6,  Os  arg£  s  tr  (Figure  3.3), 

(3.40)  A(!)  =  AV(4)  +  0(1)  . 

3.6.  The  approximation  (3.40)  for  A(£)  gives  us  the  connec¬ 
tion  formula.  It  will  be  shown  in  Section  5  that 

(3.41)  AV(Sr)  =  A£  -  2B^isin(Vw)  +  o(l)  (for  -  -  <  Y  <  ^  ) 
Therefore,  by  (3.40) 

(3.42)  AJ.=  A^  -  2B£i  sin  (Yw)  +  o(l) 

Of  course,  Af  and  A^  are  the  values  of  the  modulation  co¬ 
efficient  A(£J  at  |  =  6(e),  -5(e),  even  though  we  may  really  be 
interested  in  A(|)  for  larger  |,  e.g.,  |  =  1,  -1.  However,  in 

Chapter  6  it  will  be  shown  that,  given  a  very  slight  strengthening  cf 
assumption  (A),  Af  =  A(l)  +  o(l),  A^=  A(-l)  +  o(l),  and 

=  B(-l)  +  o(l);  thus  the  same  connection  formula  (3.42)  will  still 
hold  with  Ar=  A(l)  and  A^  =  A(-l). 

Equation  (3.42)  solves  only  half  of  the  connection  problem;  the 
other  half  is  to  show  that  Br  =  B^  +  o(l).  Given  (a"),  the 
strengthening  of  A  in  Chapter  6  Bf  =  B(l)  +  o(l)  and 
Bf  =  B(-l)  +  o(l).  Thus  it  suffices  to  show  B(l)  =  B(-l)  +  o(l).  This 
follows  from  the  well-known  WKB  formula  (2.2)  and  the  fact  that 
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B(£)  is  the  coefficient  of  the  dominant  term  of  (2. A),  for  positive 
imaginary  £. 

Mare  precisely,  we  apply  (2.2)  twice  to  £  =  i,  once  in  a 
region  R^  containing  £  =  -1  and  once  in  a  region  R^  containing 
£  =  1  (Figure  3.4) . 

The  resulting  equations  are 


z=z(i) 


q 


z=z(i) 


(A(-l)  +o(l)]e"1/£  +  [B(-l)  +o(l)]e1/e 
[B(-l)  +  o(l) ]e1/e 

[A(l)  +  o(l)]e‘1/£  +  [B(l)  +  o(l)]e1/£ 
[B(l)  +  o(l)]e1/£ 


Therefore , 


(3.43) 


B(l)  «  B(-l)  +  o(l) 
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4.  Properties  and  Error  Bounds  of  the  Kernel. 

4.1  The  conputatiorts  in  Chapter  3  depended  upon  equations 
(3. 2) - (3. 4)  and  (3.15)  for  j(r,s).  In  this  section  we  shall  derive 
those  estimates,  using  assunption  (A)  of  Section  3,  but  not  the 
restriction  (3.18)  on  6(c).  The  first  step  is  to  compute 

“i  CO 

(4.1)  j(S)e“P5  =  /  e"p%(s)ds 

4 

Assumption  (A)  suggests  that  we  approximate  ^(s)  in  (4.1)  by  y/s, 

but  in  order  to  keep  the  error  term  small  we  must  keep  the  increasing 

exponential  factor  exp(~;:s)  under  control,  and  that  requires  care 

in  choosing  the  path  of  integration  and  computing  error  bounds . 

For  the  computation  of  j(|),  the  path  of  integration  will  be  that 

shown  in  Figure  4.1  •  We  have  |  on  the  circle  between  =  -  6 

and  4r  =  5,  not  too  close  to  i.e.,  (4  j  =  6  and 

0  £  argil  £  ir  -  Q  ,  with  0  >  0  .  We  build  the  path  out  of  four 
7  o  o 

curves.  The  first  is  Aj,  a  vertical  line  drawn  from  |  =  4^  + 

A 

to  | ,  the  point  with  the  same  real  part  and  a  lower  imaginary  part, 

4^  =  mir:^!^  4i  -  —  4jj)-  This  is  chosen  so  that  Im(s),  and  hence 
|exp(-fs)|,  will  be  decreasing  along  .  The  second  curve  is  A,, 
a  horizontal  line  from  4  to  4  ,  the  point  with  the  same  imaginary 
part  and  whichever  nonegative  real  part  makes  4  lie  on  the  circle. 


-35- 


That  is,  I".  =  f\,  LiO  and  ||"|  =  6  .  This  curve  is  chosen  to  avoid 
1  1  ft 

s  =  0,  where  <p( s)  becomes  large,  without  resorting  to  an  increase 
in  |  exp(-  ps)|  .  Then  most  of  the  magnitude  of  the  integral  in  (4.1) 
will  come  from  near  £,  where  exp(-Ps)  is  biggest. 

Finally,  the  third  curve  is  the  arc  from  |  to 

=  6  exp (- it  i/4)  =  6<\/T/2  -  i&sfT/Z  ana  the  fourth  curve  is 

the  vertical  ray  straight  down  to  -i®.  These  curves  are  chosen  to 
lead  to  a  practically  fixed  point,  with  |exp(-ps)l  monotonic.  On 
A3  and  A^,  |exp(-ps)|  decreases  to  zero  exponentially  faster 
in  s  than  |  <p(s ) j  can  increase. 

Now  we  shall  approximate  <p(s)  by  y/s  and  split  the  remain  ng 
error  term  into  four  parts,  one  for  each  of  the  four  curves  Aj,  A^, 

A3,  and  A^.  By  (4.1) 

j(£)e~p^  =V/1  e"pss-1ds+2  /  e~  ps(s  <p(s)  -  Y  )s’1  ds  . 

*  i  =  1  Ai 

Now  change  variables  in  the  first  integral  and  pull  g(6)  out  of  the 
next  three  integrals 

CD  . 

(4.2)  )(£)e""p^  =  Y  /  e” u  u”  au 

PI 

+  ®(g <5))  S  /  |  e”Ps  s_1ds| 

i=  1  Ai 

+  &(1)  f  |  exp(-ps)(*(s)  -  Y/s)ds|  . 

A. 

4 
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This  holds  as  c  -*  0,  uniformly  for  r  on  the  upper  semicircle, 

U!  =  5. 

The  first  term  in  (A. 2)  is  the  incomplete  ganma  function, 

7(0, of)-  We  may  approximate  this  by  [Olver,  1974,  page  110]: 

(4.3)  r(0,  p|)=  f  e“uu_1du  ~  (p|)-1  exp(-p|)  as  e-0 

Pi 

uniformly  for  §  on  the  upper  semicircle,  not  too  close  to  i.e. 

||  |  =  6  and  0  £  arg£  *  it  -  0  (Figure  3.1  (a)). 

Once  we  show  the  other  four  terms  of  (4.2)  to  be  smaller,  this 
will  give  us  a  simple  formula  for  j(§)  • 


4.2.  The  first  of  the  four  error  terms  in  (4.2)  may  be  esti¬ 
mated  by  writing  it  in  terms  of  real  variables  and  integrating  by 
parts  to  isolate  the  contribution  of  the  path  segment  near  f , 
where  the  integrand  is  biggest. 


A 

(4.4)  J  |e-pss"1ds!  =  J  1  1  exp(-ip(£  -  t))|  |  -  it  | _1  dt 

Aj  0  1 


-1 


::  exp(-ip^1)(l  p|  r1-  \p%  -  ip(£.  -  Vi)r1exp(ip(6i  -  ?.)) 


+  ipr1 


exp(ipt)(|i  -  t)j|  -  it!-3  dt] 


The  first  term  in  the  brackets  is  e/6  .  The  exponential  in  the 


second  term  is  bounded  because  ip  <  0  and  |^  -  |^>  0-  The  deno¬ 
minator  of  the  second  tern  is  ;  p£j  . 

For  |R*  §14  we  have  f=|R  +  |i|.,  so  that  ifl  a  . 
Moreover,  for  such  |  ,  w/4  s  arg  |  s  t  -  ©q  ,  so  that 

'4.5)  S/Uf  a  sin  0q  >  0  . 

Thus,  lei  a  |  ill  sin  eo  = 5  sin  6o  . 

On  the  other  hand,  if  |  a  |  then  0<arg|  £  tt/4,  so- 

tv  i 

(4.6)  $^|||  i  COS(-ar/4)  =  l/'i'T  . 

Therefore,  ||“|  a  |  |R|  2:  |||/«/T>  j||  (sin  eQ)/2  = -|  6  sin  eQ . 

In  either  case,  the  second  term  in  the  brackets  in  (4.4)  is 

(4.7)  - 1 pf  r1exp(iP(ii - 1^)  =  iPr1!ir1&(i) 

<  (-|)(^6  sineo)"1  ©(1)  =  ©(e/5)  . 

For  the  third  term  in  the  brackets  of  (4.4),  we 

examine  the  quotient  in  the  integrand,  (|^  -  t)||  —  it  [  ~ 3  .  If 
4R*ti.  then 

||  -  tt|  *  |lm(|  -  it)|  4?.  =  ^^  . 

Therefore,  and  by  (4.5), 
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(4.8)  (ii  -  t)U  -  it  j  “  3  *  S^i"2=  =  ®<6"2) 

On  the  other  hand,  if  Li  then  we  may  use  (4.6). 


(4.9)  (et - t>i e  eteft3=  ®(eilir,9>  = 


Substituting  these  bounds  and  the  definition  of  P  =  2i/e  into  the 
last  term  of  (4.4)  shows  that 


(4.10)  IP1,-1/1  ^  exp  (ipt)(£.  -  t)\i  -  i  1 1  “  3  dt 

J0  1 

=  ©(£6_2)  f^1  ^  exp(ipt)dt  =  ©(e/5)2[l  -  exp(ip(|.  -  t^)! 
0 


=  ©(e/5)‘ 


This  completes  an  estimate  of  (4.4) 


(4.11)  f  !e“pss_1ds|  =  exp(-ip  £.)[©(  e/5)  +  ©(e/6)  +  &(e/5)2] 
'A, 


=  ©aper  exp  (-?& )) 


and  the  first  of  the  four  error  terms  in  (4.2)  is 


(4.12)  ©(g(S))  J  |  exp(- ps)s-1ds|  =  ®(g(5))(p|)  exp(-P|) 

A1 

=  o(p|)_1  exp(-  P|)  . 
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This  sunroary  error  boind  is  good  enough  for  our  purposes  because  the 
first  term  of  (4.2)  is,  by  (4.3),  approximately  y(oO  ^exp(-o'). 

It  takes  care  of  the  biggest  of  the  error  terms  in  (4.2),  because 
the  exponential  factor  in  the  integrant  is  biggest  near  |. 

4.3.  The  second  of  the  four  error  terms  in  (4.2)  can  be  inte¬ 
grated  explicitly  as  follows. 

f  |  exp(-ps)s-1ds|  =exp(-iP|i)|  /  +  1  dt j 

A2  h 

(4.13)  =  exp(-ip|'l)|lcg(t  +  |t  +  ifj  )lt= 

ft 

=  exp(-ip‘|>i) |log[(|R+  If!)  ~1(|'R  +  l|l  )]| 

To  get  a  bound  for  (4.13),  we  look  first  at  the  denominator  in  the 
logarithm.  For  |Rs  0,  f  =  |R  +  i-^lj,  s0  ^at 

t/ltl  =!&/!$  I  sm  0Q 


and 


(4.14)  |R  +  ||!  =  |f |  -  ||R|  =  (1||  +  |*Jfl(|f|2-  ||R|2) 

*  \  III"1  if/  *  III  (sin  eo)2/s 


This  bound  also  holds  if  |Ri  0,  for  then 

(4.15)  |R  +  It!  &  ltl>  it!  (sin  ©0)2/8  • 

These  two  inequalities,  together  with  (4.5)  and  (4.6),  can  give  us 
an  upper  bound  for  the  quotient  in  (4.13): 

(4.16)  (|R  +  |tl  )"1(|'a+  III )  <  (iti  (sin  e^/sf1  (2 It!  > 

=  1 1|-1  (sin 9o)"2  16  6  C  (^5  sin  ©0f  ^sin  ©0f  216S 
=  32  sin"  3  0Q  ■ 

We  also  need  a  lower  bound  for  this  quotient.  First,  if  |Rs  |^, 
then  t  =  lR  +  so  that 

III  s  III  =  6  =  |||  • 

But  |=  |R+ it^.  hence,  |R  a  |R  (Figure  4.1  (a)).  Therefore, 

(4.17)  (;R+  ltlf1(|R+  Iti  >*($«  +  ltl)'1(lR+ltl)=  1>  (sin  e//i: 
On  tie  other  hand,  if  |R  a  |i#  then 

Iti  I  =  llj-^ftl  S  max(|if3|ft)<  6  • 
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Consequently, 


III  S  |ft  +  lttl  <  2  6  • 

Therefore, 

(4.20)  (|R  +  |||  )_1(|ft+  |||  )  a  (6+  2  6f1|||  =  l/3>  (sin  0o)3/32  • 

Now  that  we  have  both  an  upper  bound  and  a  lower  bound  for  this 
quotient,  we  can  estimate  its  logarithm. 

(4.19)  |log[(|R+-  |f  |  f^f  If!)]  I 

=  max  (log  [(|R+  |||  f1(fR+  III  )L  l°g[(|R+  III  )(|R+  ill  f  *]) 

<  log[32(sin  eQf  3] 

From  (4.13),  therefore 

(4.20)  f  |  exp(-  Ps)s”*ds  |  <  exp(-ipfi)log(32  sin-3  eQ)  . 

A2 

But  is  bounded  below  |t .  For,  if  |Rs  |.,  then  by  (4.5) 

^1  “  !i  =  *  2  1 1 '  Sin  Go  =  2  6  sin  eo  ‘ 


If  then  by  (4.6), 

Therefore,  (4.20)  gives  us  room  to  keen  the  integral  smaller 
than  an  incomplete  gamma  function: 

(4.21)  f  |exp(-ps)s_1ds|  =  ${exp(-ip£\)) 

A  1 

2 

=  6(exp(ip(y6  sin  eQ)))exp(-ipei) 

*  S(exp(-6/e))|exp(-p|)i  =  o(p£ )_1exp(-  p|)  . 

This  shows  that  the  second  of  the  four  error  terms  of  (4.2)  is  smaller 
than  the  first  term  of  (4.2). 

4.4.  The  third  of  the  error  terms  in  (4.2)  if  readily  estimated 
by  repeating  the  computation  of  (4.21): 

(4.22)  J  |  exp(-  ps)s_1dsi  *  I  exp(-  p|)|  5_1  f  |ds| 

A3  A3 

*  expt-ip^JS'^irS/^  =  exp(-  ip£.) 

=  otp^^exp^pg)  . 
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The  last  of  the  error  terms  in  (4.2)  is  the  reason  for 
our  assumption  that  <?{£)  =  ©(ec  ^  1  )  as  4  — -i®  (Chapter  2)  . 
The  other  end  of  A  is  4b  =  5  exp  (-  tr  i/4);  there,  <?(£,)=  G- (1/6 ) 
by  assumption  (A).  Accordingly, 

lA  |  exp(- ps)(v(s)  -  Y/s)ds!  =  ©(1/5 )J  |exp(K|si  -ps)ds| 

4  A4 

00 

=  ©(1/5)  /  exp((<  +  ip)t)dt  =  (1/5 )(-<  +  2/t)-1exp((*:  -  2/e  )  6) 

6 

=  ©(e/6)  exp (-2 6/e  ). 


But  the  4  which  we  are  considering  here  have  positive  ima¬ 
ginary  parts,  so  that  p£  has  negative  (or  zero)  real  part,  so 

exp(-26/e)  =  ®(exp(-  p  | ))  exp(-  26/e  )  =  o(exp(- p|  ))  . 


Therefore, 

(4.23)  f  |exp(-ps)Ms)  -  Y/s)ds|  =  (P4)-1exp(- p|)©(exp(- 26/e  )) 

A4 

=  o(P4)'1exp(-p4)  . 


Now  substitute  into  equation  (4.2)  the  equation  for  its  dominant 


term,  (4.3),  and  the  equations  for  its  four  error  terms,  (4.12), 
(4.21),  (4.22),  and  (4.23).  The  result  is: 


^  1  ®  1  . 

(4.24)  j(4)e"P^  =  f  e'pS(V/s)ds  +  S(g(  S)(P4 )  e“P^ 

I 

+  afexp^fi/eHfpfcfV^ 

—i«o  - 

=  (1  +  &(h(e  )))  f  e"ps(V/s)ds 

t 

=  (Y  +  o(l))(P4f1e"p^  . 

This  holds  uniformly  in  4  on  the  arc  |  4  ]  =  6 ,  0  <  arg4  s  ir  - 

(Figure  3.1  (a)). 

Substituting  this  into  (2.12)  gives  an  equation  for  j(4,  s) 

(4.25)  j(4  ,  s)e~ps  =  (1  +  <&(h(e))jY(4,s)e'ps 

=  (V  +  o(l»[(ps)"1e"ps  -  (P4)"1e“p^]  . 

This  holds  uniformly  in  4  and  s  on  the  arc  U  !  =  I  sl  = 
with  0  *  arg  %  <  it  -  eQ,  0  <  arg  s  <  tv  -  0q  (Figure  3.1  (a)). 

4.5.  We  still  lack  an  estimate  for  j(4,s)  with  4  or  s  near 
4^,  where  (4.3)  fails  because  4^  is  a  sector  boundary  of  that 
approximation.  By  (2.12), 

e 

(4*26)  Htr.^)sxp(-P4^)  =  fT  e‘Ps<p(s)ds  . 
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and  by  (A) , 


(4-27)  *<s)=  (Y/s)  +  (g(s)/s)  . 

with  g(s)  ->  0  as  |s|  ■*  0. 

In  the  lower  half  s-plane,  ps  has  positive  real  part,  and  the 
exponential  is  small,  but  has  a  possible  branch  cut  along  the 
negative  imaginary  axis  (Figure  2.1).  Write  the  upper  semicircle 
(Figure  4.2  (a))  as  the  sum  of  curves  L  and  F  (Figure  4.2  (b)). 
Ihen 

(4-28)  J(£r*£  *)  exp  =  Y  ^e'^s^ds 

+  J  e  pSip(s)ds  +  f  e  Psg(s)s  *ds  • 
r  L 

The  integrand  of  the  first  integral  in  (4.28)  has  only  a  simple  pole, 
(4.29)  Y /  e“PS  s_1ds  =  -2:riY  . 

V 

We  shall  show  that  this  is  the  dominant  term  of  (4.28);  the  second 
term  has  a  decaying  exponential  and  the  third  term  includes  g(s) , 
which  is  small. 

The  second  term  of  (4.28)  would  normally  be  shown  small  by 
Jordan's  lemma,  but  the  standard  forms  of  Jordan's  lemma  do  not  cover 
this  case.  Break  r  into  three  parts,  as  shown  in  Figure  4.2  (e) . 
Because  •f  is  holoirorphic ,  assunption  (A)  implies  that  on  F^, 
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(4.30) 


d<?(s)/ds  =  G(s-2)  =  &(5-2)  • 


Therefore  we  can  get  a  useful  error  bound  by  integrating  by  parts. 


-1  -ps  .  .  i’a  -1  -Ps 

^  in  *  I  c'  \  :  _  /n 


f  e-ps  <p(s)ds  =  -p-ie  ?( s)L“  -  p  Ae  ra?(s)  , 

r,  ^b 


(4.31) 


+  p-1  J*  e-ps  <?'(s)ds 


-  ©(e/6)  . 


For  the  second  part. 


Im(s)  s  -  5//2 


i 


so  that 


(4.32)  f  e~p%(s)ds  =  ®<exp(-«/Z(6/e ))) 
r2 


/  |  <p(s)ds|  =  o(e/S)  • 

r2 


For  the  third  part,  we  use  the  assumption  that  <p  grows  at  most 
exponentially  as  s  —  -  ice,  and  assumption  (A)  for  s  —  0  . 


(4.33)  |  f  e"ps<p(s)ds[  s  25 ~i  f  exp ((/c  -  2/e)t)dt  =  ©(e/6) 


-1 


Combining  (4.31),  (4.32),  and  (4.33)  gives  a  bo  md  for  the  second 
term  of  (4.28) 


m  : 
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im 


i 


(4.34) 


f  e~^S  <p( s)ds  =  &(£/6)  . 


In  the  third  term  of  (4.28),  the  exponential  factor  in  the  integrand 

I 

can  grow.  We  shrink  the  circle  in  L  from  a  radius  i  to  a 
radius  e  (Figure  4.2  (b)  and  (d));  then  on  the  resulting  new  curve 
L,  the  exponential  satisfies 


(4.35) 


exp(-  ps )  |  s  exp  (-pie)  =  e  . 


Thus  the  third  term  of  (4.28)  is,  where  L,  is  the  circular  part  of 
L  and  L 2,  the  straight  part. 


(4.36)  f  e"’pSg(s)s‘"1ds  =  f  e~^sg(s)s  *ds 
L'  L 


_  r  e-psg(s)s  *ds  +  f  e  °sg(s)s  ^ds 
^L1  L2 

=  G(g(£))+  /  e_Ps  g(s)  s-1ds  . 

L2 


In  I2,  the  last  integral  of  (4.36),  g(s)  is  small  where 
exp(-ps)  is  biggest.  However,  standard  forms  of  Watson's  Lemre  do 
not  apply. 

Let  g^(s)  represent  the  value  of  g(s)  on  the  right  edge  of 
the  negative  imaginary  axis  and  g_(s)  the  value  on  the  left.  Thai, 
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-ioo 

(4.37)  I  =  f  exp(-  ps)(g  (s)  -  ,  (s))s_1  ds  . 

-i  e 

With  t  =  is  -  •  and  p(t)  =  g  (s)  -  g^(s) , 

(4.38)  I  =  e-2  /  exo(-2t/£  )(t  r  £)"1p(t)dt  . 

0 

Now  we  break  this  integral  into  three  parts.  First,  consider 
til-  Then  because  <3(3)  =  <t(e:<p(/c  ,  si  ))  for  s  —  -  i®, 
also  p(t)  =  ©  (exp{,<  t)) .  Therefore, 

<®  .  03 

/  exp(-2t  /e)(t  +  e)  p(t)dt  s  f  exp((/c  -  2/e  )t)dt 

1  1 

(4.39) 

=  (2/e  -  k  )-1  exp(K  -  2/e)  =  &(e  )  exp  (-2/e  )  . 
Second,  consider  ^6  s  t  s  1  ;  in  this  interval, 

I  p(t) I  s  2g(jt|  +  E)  s  2g(l  +  e  )  <  2g(2)  , 

(t  +  e  f 1  &  (|o  +  e  f1  <  26-1  , 


(4.40) 

=  ©(e/6)  exp  (-  26/e  )  . 


r  -1  1  1 

/  exp(-2t/e)(t  +  e)  p(t)cts  ©(6)  f  exp(-  2t/e  )dt 

‘Is 
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Lastly,  consider  0  <  1 1  -^6  .  Here, 

|p(t)|  <  2g( !  1 1  +  £)  £  2g(j5  +  e)  =  &(g(S))  » 

,  76 

exo  (-  2t/s  )(t  +  £)“  p(t)dts  S(g(6))£_I  f“  exp(-2t/e)dt 

"0  0 

(4-41) 

=  s(g(o))  • 

Finally  coiTbine  (4.39),  (4.40),  and  (4.41)  to  bound  ^ 

(4.42)  I2  =  S(g(6))  +  ®(e/5)exp(- 2  6/e)  =  o(l)  . 

From  (4.42),  (4.36),  (4.34),  (4.29)  and  (4.28), 

(4.43)  j(|r,|£)exp(-p|.)  =  -  2-rr iV  +  ®(g(6))  +  fe(e/5) 

=  -  2ir  iV  +  o(l)  . 

4.6.  Besides  computing  we  must  find  out  how 

close  it  is  to  By  (4.28)  and  (3.14), 

(4.44)  i\r,Z£))exp(-pZ£)  =  j^e"ps g(s)s"1ds 


+  J  ,  e”ps  g(s)s  ^ds  . 


By  (4.36)  and  (4.42),  the  second  term  is 

(4.45)  f  e_psg(s)s_1ds  =  ©(g(e))  +  ®(g(6))  +  G(e/6 )  exp  (-  25/e  ) 

L 

=  $(h(e))  • 

For  the  first  term  of  (4.44),  we  use  equations  (4.31),  (4.32),  and 
(4.33)  with  g(s)s~*  in  place  of  <p(s).  Because  g  is  holomor- 
phic,  on  r1# 


(4.46)  dg(s)/ds  =  &(g(5))/S 

®(g(6)6'2)  • 

Substitute  this  into  (4.31)  to  find  that 

(4.47)  /  e~Psg(s)s_1ds  =  ©(e5_1g(S))  . 

rl 

If  g  is  substituted  for  <p  in  (4.32)  and  (4.33),  we  get 

(4.48)  /  e"psg(s)s_1ds  =  o(eS-1g(5)) 

r2 

(4.49)  f  e~Psg(s)  s~*ds  =  ©(£  5“*g(5))  . 

r3 

Therefore, 
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(4.50)  (JCe^  e£)  -  iV($r.£je))exp(-F$2)  =  S(h(e))  +  &(£S_1g(S)) 

=  &(h<£))  =  fe(h(£))jY(|r,ei)exp(-peA.  ))  . 

4.7.  Equation  (4.43)  leads  to  a  more  general  equation.  Add 
j(C.£r)  and  use  (2.12)  and  (4.25), 

(4.51)  j(44x)exp(-pei)=  (Y  +  oU»[(p4rf1expf-p£r) 

-  (PSf^expf-pl)  -  Zir i] 

as  e— 0,  uniformly  in  4  on  the  arc  |||  =  5,  Os  arg|  s  it  -6o 
(Figure  3.1  (a)). 

This  has  still  not  dealt  with  j(4,s)  for  £  and  s 
both  near  4  .  Next  we  shall  find  a  bound  for  j(4,4?) 
for  |  near  % ^  ;  i.e.,  ir  -  0Q  s  arg  4  s  it  .  By  assumption  (A), 
<p(s)  =  ©  s-*)  uniformly.  Therefore,  if  we  perform  the  integration 
(2.12)  on  an  arc  of  the  circle  |  s|  =  6  and  set  t  =  Im(s)  (Figure  4. 3), 

(4.52)  j(4,4^  )exp(-  p4^)  s  |  exp(-  ps)<?  (s)dsj 

& £ 

=  &(1)  J  I  exp(-  ps)s_1ds| 

h 

I  .1 

=  &(1)  /  1  exp(-ipt)(52  -  t2)  2  dt  • 

0 
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2  2 

Now,  because  -  Gq  s  arg  s  <  it  ,  we  can  bound  (6  -  t  ) 

as  follows: 

Ost  =  Im(s)  =  6  sir.  arg  s  =  6  sin  (tt  -  arg  s)  s  6  sin  6q 

(4.53)  62  £  (62  -  t2)  i  5Z(1  -  sin2©o)  • 

Substitute  this  into  (4.52). 

(4.54)  j(| ,  4  „)exp  (-  ,)  =  ©(6-1)  f  1  ex p(-ipt)dt 

*  1  o 

=  &(e/6)exp(-ip|.)  =  <^P4  f1  exp  (- P4  ) 

But  by  (4.51)  this  bound  also  applies  to  any  |  on  the  arc 
l4l  =  6,  eQ  £  arg 4  s  it  -  0Q  (Figure  4.4)  •  Consequently, 

(4.55)  j(4, 4^)exp(- p4^)  =  &(P4)_1exp(- F4)  as  e  -  0 

uniformly  for  4  on  the  a.c  |||  =  5,  6Q  £  arg 4  s  it  (Figure  4.4). 

Now  use  this  equation  twice,  to  subtract  j (s , 4„)exn(-p4f ) 

from  j(:„  4^  )  exp  (-  P4^  )  .  By  (2.12),  for  4  and  s  as  in  (4.55), 

(4.56)  j(4,  s)exp  (- ps)  =  (5(P4)-1exp  (- P|)  +  e(ps)-1exp  (- ps)  . 
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These  are  both  increasing  exponentials  in  (4.56)  .  If  we  restrict 
£  and  s  to  the  left  half  of  the  semicricle  with 
it  s  arg|  £  arg  ss  it  (Figure  3.1  (c)),  then  the  first  term  dominates, 
so  (4.56)  reduces  to 

(4.57)  j(£,  s)exp  (-  ps)  =  ©(p|)-1exp  (-  p£)  as  e—  0. 

4.8.  To  complete  the  proof  of  the  estimates  (3.17),  (3.18) 
and  (3.19),  which  were  essential  to  Section  3,  observe  that  (3.17) 
is  the  same  as  (4.25).  For  equation  (3.18),  add  (4.25)  with  s  =  f 
to  (4 . 50)  . 

To  get  equation  (3.19),  we  must  modify  the  conputations  in 
section  4.7.  First,  subtract  (2.12)  from  itself  as  follows. 

(4.59)  [j(|4jg)  ~  1  Y(!*!jg)]  exp  (-  =  f  exp  (-  ps)g  (s)s_1ds 

* i 

£  &(g(5))  f  |  exp  (-  ps)s_1ds| 

The  last  integral  also  appears  in  (4.52,  and  was  bounded  in  (4.54), 

by 

(4.60)  ©(pSJ^ex p(-p|)=  © (£  6-1)  exp  (-  P£ )  . 


There  fo-e, 


(4-61)  [](£,£*)  -  Jy(£.S£)]exp(-F€£)  =  &(S)  E  5-1exp(-p£)  • 


Again  subtract  (2.12)  from  itself,  and  apply  (4.61)  to  the  result. 

[Hi,  s)  -  jY(|,s)]  exp  (-  ps) 

(4.62)  =  [Hi,i*)  _  ^  (£»4^  ^  3XP  (”  )  -  [Hs,4^)  ~j  i2,i^ )]  Qxp(~  »  ) 

=  ©(g)  e  6-1  exp  (-  p£ )  +  ©(g)  e  6  1  exp  (-  ps) 

But  equation  (3.19)  is  required  to  hold  only  for  sarg£  *  arg  s  s  tt 
(Figure  3.1  c);  there,  exp(-  ps)  s  exp  (-  ).  Thus,  (4.62)  reduces  to 

(4.63)  [j(£.  s)  -  Jv($,  s)]exp(-ps)  =  ©(g)  e  5-1exp(-P£) 
which  is  equation  (3.19)  • 

4.9.  In  Section  6,  where  connections  are  made  for  t  on  the  real 
line,  we  will  need  a  bound  on  j(l,s)  for  Z  and  s  real,  with  the 
same  sign  Since  <*>(s)  is  holomorphic  and  q>( s)  =  ©(s-^)  as  |  s|  —  0, 

for  all  s  not  on  the  negative  imaginary  axis,  the  bound  d<j(s)/ds 

=  ©(s  j  holds  for  real  s  —  0.  Therefore  we  can  integrate  (2.12) 
by  parts: 
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,,s)e‘Ps  =  -P-1  e'’pt*(t)|?  +  p'1  / e_pt  ?'(t)dt 


-1  -2 

=  $(et  l)  I  +  S(z)  J  t  dt 


-1  .6 

=  &<et  )  is 


If  6,  s  ||| ,  I  s|  s  5i  +  j,  then  this  implies  that 


(4.65) 


id.  s)  =  &(eS7  ) 
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5 .  Completion  of  Close  Connection. 

5.1.  We  shall  compute  c(r)  when  q(z)  =  z'log'z  (■  '■  -1) 
near  z  =  0  to  show  that  assumption  (A)  of  Section  3  holds  for 
such  q(z) .  In  this  subsection,  we  use  the  notation  log'  z  =  (log  r.) 
and  asymptotic  relations  hold  as  z  -*■  0;  i.e.,  f  -*•  0. 

(5.1)  q'  (z)  =  zv_1  log4  z  (v  +  a  log"1  z) 

*  z  *q(z)(v  +  p.  log'1  z)  =  z  ^(zHv  +  o(l)  )  . 

Substitute  this  into  the  definition  (2.7)  of  c(0. 

(5.2)  ^=jq2q'  =jz1q1(v  +  ji  log  1  z) 

=  ~  (z  q)"1  (v  +  o(l) )  . 

To  compare  this  with  assumption  (A)  we  must  write  zq  in  terms 
of  4  .  Integrate  the  definition  (2.3)  of  L  by  narts  and  use 

(5.1), 

.2 

(5.3)  |  =  z  q  -  f  q'  z  dz  = 

0 

-1  rZ 

=  z  q  -  v  4  -  8  (log  z)  J  |  q(z')dz'  | 

0 
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To  see  that  the  last  integral  In  (5.3)  is  of  the  same  order  as 
set  z  =  r  exp  (ia)  with  0  <  r  <  1  and  write  |  4 !  as  : 

(5.4)  !  4  I  =  I  f  sV  (log  s  +  ia)^  ds  =  1  J  sv  log^  s  ds  j 

0  0 

+  ©  (log  *  r )  f  |  sv  log^  s  ds ! 

0 

=  (1  +  ©  (log  1  r) )  f  |  s v  log  ^  s  ds  | 

0 

i 

=  (1  +  &  (log  1  r) )  J  |q(z')dz'| 

0 

Therefo  e,  (5.3)  says  that 

> 

(5.5)  zq=|(l  +  v+  ©(log-1  z))=£(l+v  +  o(l))  . 

Now  (5.2)  and  (5.5)  imply  assumption  (A): 

(5.6)  *>=  j(l  +  v)"1  v  e_1(l  +  ©(log-1  z)) 

*  x  (1  +  vf1  V  4'1  (l  +  0(1))  . 

With  (5.5)  and  the  definition  of  g(z) ,  this  also  gives  us  g( 

log  4  +  o(l)  =  log  z  +  log  q  =  (1  +  v  +  o(l) )  log  z 

(5.7)  g(6)  =  ©  (log-1  z)  =  ©  (log-1  i)  =  ©  (log-1  5 )  , 
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5.2.  We  need  the  connection  formulae  for  A', 


that  is ,  for 


(5.8) 


q(z)  =  z 


(v  >  -1) , 


we  must  obtain  the  equations 


(5.9) 


(a)  AY(1)  =  AY{6)  +  o(l) 

(b)  AY(-1)  =  AY  (-6)  +  o(l) 

(c)  AY  (1)  =  AY(-1)  -  BY(-1)  2i  sin  (Ytt)  +  o(l) 


where 


(5.10) 


Y  =  j  v  (v  +  If1 


(-co  <  Y  <  y) 


When  (5.8)  holds,  the  solution  w(z)  to  (2.1)  is 
explicitly,  in  terms  of  the  Hankel  functions  and 


(5.11) 


w=  (4/e)  +  C2  4XH^2)  (4/e) 


where  Cj  and  C2  are  arbitrary  constants  and 


(5.12) 


^=j(v  +  l)1=  j  -  Y  . 


•i 

i 


The  connection  formulae  will  follow  from  (5.11),  (2.2)  and  the 
appropriate  formulae  for  the  Hankel  functions.  It  is  known  [Olver, 
1974,  pp.  238-2391  that 


(5.13)  (a)  H^(4/e)  ~  e*^2  %  exp(i4/e) 

for  -7r  +  eQ  <  arg  %  *  rr  ,  as  |  4/e  |  >  (6/e)~cc 

(b)  H^2)  (  4/e)  ~  A?  e1/2  4_3/2  exp(-i4/e) 

for  -rr  s  arg  4  s  tt  -  eQ  ,  as  |  4 /el  a  (5/e)  -  » 

(o)  (4/e)  =  (4  exp  (-rr i)/e)  exp  [\  rri)  + 

+  (4  exp  (~~i)/e)  2  cos  (X~) 


with 

=  (2/rr)1,/2  exp  (  -  in  (  \  +  1) ) 

A2  =  (2 /rr )1/2  exp  ( ^i77  (  \  +  |- ) ) 

We  can  substitute  (5.13)  (a)  and  (b)  into  (5.11),  and  this 
be  compared  with  (2.2)  because 

(5-l5)  q  l//2  =  (constant)  4X‘~l//2  . 

by  (5.8),  and  (5.9)  (a)  follows. 


can 
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Equation  (5.13)  (b)  does  not  cover  the  negative  real  nurrbers 
which  (5.9)  (b)  and  (c)  need,  but  we  can  extend  it  by  using 
(5.13)  (c). 

(5.16)  (|/e)  ~  2  i  ^  £  1/72  cos  (Xtt)  exp  (i|/e)  + 

+  iAj  e1//2  4  1,/2  exp  (\n i)  exp  (-i4/e)  . 

Now  equation  (5.9)  (b)  follows  from  (5.16),  (5.13)  (a)  and  (5.11)  . 
For  the  last  part  of  (5.9)  we  combine  (5.16),  (5.13),  and  (5.11) 
once  more. 

A  (1 )  =  CjAj  e1/2 

(5.17)  A  (— 1)  =  ClAl  e1/2  +  C2  ZiP^&/Z  cos  (XT r) 

1/2 

B (-1)  -  C^iAj  e  exp  (XttI) 

Therefore,  by  (5.17),  (5.14),  and  (5.12), 

A (1)  =  A(-l)  -  C2  ZIA^  £^2  cos  (X'7T) 

=  A(-l)  -  B  (-1)  (Aj  s  /  exp  (\tt1)  )  *  2  e  ^  cos  (\7T ) 

=  A(-l)  -  B(-l)  21  cos  (X- ) 

3  A(-l)  -  B  (-1)  21  sin  (YTT) 
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5.3.  In  Section  3.2,  we  had  claimed  that  the  method  there  use., 
could  be  carried  out  more  precisely  to  compute  A()  instead  o! 
just  bounding  it.  To  do  that,  take  r ^  =  4.,  and  j  -  _  arg  •'  _  - 

lit  =  5  (Figure  3.1  (d))  and  define 

r€ 

I  a  (  4)  =  J  a  (s)  j  (| ,  s)  exp  (  p  |  -  p  s)  <?(s)  ds  . 

h 

In  Section  3.2,  we  had  shown  that 

£ 

(5.18)  ||I  a||  *  $(e/6)  . 

00 

With 

(5.19)  a,  (4)  a  a  (4 )  -  a  exp  (p|  -  p4  .)  .  . 

1  l  l 

equation  (3. 7)  may  be  written  as 

(5.20)  ax(4)  =  J  (4  ,  4i)  exp  (  P4  -  P4^)  +  I£  a  (|) 

To  use  (5.18)  to  contract  (5.20),  we  need  Ira^(I)  and,  in  turn, 
Irexp(pi).  By  (3.4), 
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By 


I  exp  ( p4)  =  I  exp  (  ps)  j  (4  ,  s)  exp  (  p|  -  ps)  ds  = 

£  i 

i  A 

=  ®(p4)  /  I  exp  (  ps)  <p(  s)  ds|  . 

(A),  <?(s)  =  S(6-1)  .  Therefore, 


-1  r s 

(5.21)  I  exp  (P4)  =  ©(6  P4 )  J  !  exp  (  ps)  ds  |  . 

e  4 


Set  s=  6  exp  (l("-9))  (Figure  5.1)  and  use  the  fact  that 
sin  «  >  j  for  '  to  compute 


,  X  exp  ( p4)  =  ®(e5_1)  f 


(i-rr  -arg  4) 


exp  (-2  6  e  *  sin  0)  de 


(5.22) 


i  i 

c-lx  r  _  ,  r  -1 


<©(e5 ~l)  f  exp  (-6  e_I  0)  de  =  &(e/6)  . 

0 


Therefore,  if  we  apply  I  to  (5.19)  we  will  get 


(5.23)  I£  a  (4)  =  I£  ax(4)  +  exp  (-p4^)  Ic  exp  (p4) 

=  I£  a:(4)  +  ®(e/6)2  . 

If  we  substitute  this,  (5.18),  and  (3.4)  into  (5.20),  we  pet  a 
bound  on  a^('). 
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(5.24) 


aj  (!)  =  Me/6)  +  [Me/'6)2  +  Ie  a2  (i)] 

=  &(e/6)  +  Me/6)  ||at  ;| 1  . 

i  00 

Therefore , 

(5.25)  ajd)  =  Me/6)  . 

This  holds  uniformly  for  4  in  the  left  quarter- circle 

|-7r  <  arg  4  -  77  (Figure  3.1  (b))  .  Now  substitute  (5.25)  back 

into  (5.20),  using  (5.18)  again. 

(5.26)  ajd)  =  Bi  Ut.i£)  exp  (P4-  pZ£)  +  ©(e/6)2  . 

If  we  require  -  ar9  £  -  77  “  ©q  (Figure  3.1  (d) )  ,  then  we  can 

use  equation  (3.3)  for  j(|,4^)  . 

(5.27)  a1(g)  =  Bj(Y  +  o(l))  [(pi^"1  exp  (  p|  -  P4f )  -(plf1  - 

-  27t  i  exp  (p4)]  +  Me/6)2  . 

But  for  s  arg  %  s  tt  -  9  0  , 

exp  (p|)  =  ©  (e/5)2  , 

so  the  second  term  in  the  brackets  dominates  and  we  are  left  with 
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(5.28) 


aj(|)  =  -B  (Y  +  o(D)  (pi)  • 

Moreover,  the  first  term  in  (5.19)  dominates  for  exactly  the  same 
reason.  This  gives  us  the  equation  we  had  sought: 

(5.29)  a(4)  =  -B/(Y  +  0(D)  (pi)"1  as  e-  0 

uniformly  for  j77  -  arg  4  :£  tt  -  ©Q  ,  |  4  |  =  6  (Figure  3.1  (d) )  . 

In  terms  of  the  original  modulation  coefficient  A  (4)  ,  (5.2  5)  and 

(5.29)  say  that 

(5.30)  A(4)  -  At=  ©  (e6_1)  exp  (-pi)  = 

=  &(  pi)  1  exp  (-pi) 

uniformly  for  i  in  the  upper-left  quarter- circle  ,  ]  4 1  =  6  , 

jtr  <  arg  |  £  7T  (Figure  3.1  (c) )  .  Moreover, 

(5.31)  A(S)«  -B  (Y  +  0(1))  (pi)-1  exp  (-pi) 

uniformly  for  i  in  part  of  the  same  quarter-circle,  the  arc 
III  =  6,  i-7T  £  arg  i  S  n  -  ©0  (Figure  3.1  (d) )  . 

Sine*,  exp  (-pi)  increases  as  Im  (i)  increases,  this  says  that 
A (4)  intreases  in  magnitude  as  i  moves  from  4  around  to 

X/ 

4  =  16  ,  the  halfway  point  between  4  and  4  (Figure  3.1  (c) )  . 

&  r 
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For  £  bounded  away  from  £  by  a  fixed  angle ,  A(4)  increases 

1/ 

almost  exponentially  in  (6/e)  as  £-*-0  and,  hence,  (6/e)— »  • 

However,  our  choice  of  5=  6(e)  may  typically  have 

(6/e)~  log  |  log  ej  ,  so  this  "exponential"  increase  in  A  (4)  may 

actually  be  rather  slow  because  4  approaches  zero  so  fast. 

For  4  very  near  to  4,.  equation  (5.31)  does  not  apply ,  but 

l 

we  can  use  (5.30).  If  4  stay  near  enough  to  4,  ‘  for  example  , 
if  the  angle  between  4  and  £  be  no  more  than  (e/6)  ,  then  the 

Jj 

right-hand  side  of  (5.30)  is  small,  so  that  A(4)  is  approximately 

equal  to  A.  . 

At 


-Ik- 


6.  Connection  from  Near  Zero  to  Far  Away. 


6.1.  We  have  obtained  formulae  connecting  the  modulation 


coefficients  A(0  and  B(\)  at 
£  =  £  =6.  But  6  =  5(c)  -*  0  as 


=  -5  to  those  at 


0,  and  we  might  really 


want  to  connect  points  which  are  much  farther  from  zero,  e.g. , 


£  =  Ct  =  -1  and  5  =  f,_  =  +1  or  even 
£  =  CR  =  exp(c_1) . 


=  -106  c  ^  and 


The  obvious  approach  to  connecting  the  coefficients  at  £^ 
to  those  at  r  is  to  connect  £,  to  £ ,,  along  the  negative  real 
axis,  C,  to  f  by  the  methods  of  Sections  2  through  6,  and  then 
Cr  to  £R  along  the  positive  real  axis.  The  WKB  solution,  A(0  ' 
constant,  B (•" )  -  constant,  is  well  known  to  be  valid  along  the 
negative  real  axis  and  along  the  positive  real  axis  (subject  to  some 
integrability  and  boundedness  conditions  at  +  “)  so  long  as  £  is 

bounded  away  from  the  possible  singular  point  £  =  0.  This  is  good 

enough  to  connect  from  points  like  =  exp(e  to  fixed  points 
like  =  1,  but  we  still  need  to  push  on  to  £r  =  6(c) . 

This  raises  the  question,  how  close  to  zero  is  the  WKB  solution 
really  valid?  If  we  could  extend  it,  for  real  nunters,  close  enough 

to  zero,  then  fT  and  £„  would  be  connected  to  £.  and  £  ,  and 

the  problem  would  be  solved.  We  shall  now  show  how  to  carry  out 
such  an  extension  in  most  cases. 

One  approach  is  to  imitate  the  method  of  Langer  (1931) . 

Although  he  studied  only  the  case  where  q(z)  ~  zv  at  zero,  and 


ir 


obtained  an  asymptotic  formula  good  at  z  =  0,  his  method  can  also 
be  applied  to  the  general  case  to  get  a  weaker  result.  For  q(z) 
as  was  described  in  Section  2,  and  for  real  '  in  R. ,  it  yields 

(6.1)  lAU)-A(-l)  !  =  ©(e)  supt  j  (di/dr)-2i2|  :  j  $  |  <£)  } 

|B(C)~ B(—  i)  |  =  ©(e  )  sup  {  j  (d$  /d=,  ) -2;2  j  :  |  s,  j  >_i  (e  )  } 

uniformly  for  £  <_  -  a  (e) 

(6.2)  i  A  (C) -A  (1)  |  =  ©(:)sup{  |  (d?/dU-2;2|  :  |  ;|>^  (e)  } 

|B(5)-B(1)|  =  &(e)sup{  I  (d^/dc)-2E2|  :  U,  >A  (c)  } 

uniformly  for  ?  (z)  . 

Thus,  Langer's  method  can  extend  the  WKB  solution  to  a  distance 
A(c)  from  zero,  with  A(e)  ->■  0  as  c  ■*  0,  Drovided  that  the 
right-hand  sides  of  (6.1),  (6.2)  approach  zero  as  e  -*•  0. 

How  useful  is  this?  Suppose  that  assumption  (A)  of  Section  3 

holds  in  addition  to  the  assumptions  of  Section  2.  Then  because 

-1  *2 

t  =  ©(5  )  and  «J>  is  holomorphic,  (d$/d£)  =  ©U  )on  this  real 

line;  therefore, 

(6.3)  (d4>  /d£ )  -  2<p  =  ©(?  )  as  |  S  |  -*■  0  ,  uniformly. 


Thus  the  condition  for  this  extension  method  to  work  is 


(6.4) 


e  A  (e ) 


0 


as 


£  -►  0  . 


-2 


ct  1 

This  allows  us  to  take  A(t)  =  e  with  :i<~2'  ror  example. 

But  such  '  ( ■  )  will  not  be  good  enough;  we  must  extend  the 
\1KB  solution  still  closer  to  zero.  As  noted  in  Section  3.3,  a  typical 
interesting  ;(‘)  would  have  5/e  =  ^  log | log  '!;  while  if 
A(0  =  d  then  '  =  (25)  ' (logj log  5|)-cx  5. 

6.2.  To  connect  from  -  =  -1  to  r.  and  ■"  to  =  1  we 

i.  r 

will  again  use  the  integral  equation  (2.13).  Besides  the  assumptions 
of  Section  2,  we  will  use  a  weaker  form  of  assumption  (A),  Section  3 

(A')  }  (C)  =  ©(C*1)  as  |c|  0,  uniformly. 

Since  <p  is  holomorphic,  it  follows  that  for  real  £  , 

(6.5)  <$  1  (C )  =  ©<5~2)  as  |d  -*■  0. 

We  shall  start  by  connecting  between  ±5  1  (e  )  and  ±5  2  (e  )  , 
where  e  <_  5  ^  <.  5  2'  and 

-1  -12 

(6.6)  t  5  ^  (log  (5  25  ^  ))  ->-0  as  e  -*■  0  . 

Equation  (6.6)  keeps  5^  and  5^  from  being  too  far  a- 
part  for  us  to  control  integrals  ranging  between  5  ^ 
and  5  2* 
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Write  equation  (2.13)  and  its  equivalent  for  posi- 


tive  5 

,  equation 

(3.5)  , 

as  follows: 

(6.7) 

A(c)  -  A± 

=  Bx  j 

(^  ,C1)exp  (-p^1)  +  TfA(;) 

where 

(6.8) 

T  A  (5) 

£ 

=  f 
?1 

<p(s)  A  (s)  j  (5  ,s)  ds 

and  whe 

re 

(6.9) 

> 

i-* 

n 

> 

H* 

'  B1  " 

8(5^  '  m 

When 

«1  =  +  51 

,  we  take  ^  <  5  <  5  2  ; 

when 

51  *  "  51 

,  then 

-  <$2  -  5  -  ci  • 

We  can  bound  by  means  of  (A'),  (6.5),  and  (4.65).  Let 

represent  the  sup  norm  on  the  interval  over  which  r.  ranges. 


lv(5)|  i  II  AH-  115  II.  IlMlI.  I  |s_1as| 

51 

(6.10) 

<  II  All.  ®(e/61)  log(62/61)  =  j|  A  ]  |  ^  o  (1) 

Therefore  T^  is  a  contraction  for  sufficiently  small 
e,  and  (6.7)  has  a  bounded  solution.  Substitute  (6.10) 
and  (4.65)  into  (6.7) 

(6.11)  AU)  -  Ax  =  ®(e/61)  +  ||  A||  _  ©  (e  5~  Xlog  (<5  2 /6  )  ) 
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Therefore  Jj  Ajj  is  bounded  independently  of  z,  and  hence, 
(6.12)  h(i)  =  h1  +  ©Up"1  log  (5  2  /  S  )  =  A1  +  o  (1) 


uniformly  for  •;  between  5^  =  ±  5^  and  ±  52" 

To  conpute  the  other  modulation  coefficient,  B (!) ,  integrate 
equation  (2.8)  by  parts  and  use  equations  (A'),  (6.12),  (6.5),  and 

(6.6)  and  the  fact  that  the  integral  is  corrputed  on  a  real  interval . 


(6.13)  B(£)  -  B.  =  /  ,Ms)A(s)epS  ds  = 

«1 


=  Ax  p"1  epS  *  (s) 


£ 

c 

1 


/  ♦'  (s)eps 

«1 


ds 


+ 


£  DS 

+  /  *  (s)  (A  (s)  -  A1)eps>  ds 

C1 

5  -2 

=  ©(e/6,)  +  ©(e)  /  s  ds  + 

-1  ^  - 1 

+  ©U5,A)  log  (6  2  /«  . )  /  s  xds 

«1 

=  ©U /6  ^)  +  ©U/5^)  +  ©(e5^)  (log  (6  ^ /6  ^)  ) 
=  o(l) 
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;  fa  (  =  )  =  :-4  *  o  (i) 

“  ■  • .  0  IT  ~  O  ‘il  a.  1 ."}  L  O  o  JT  O  V  O  ^  Ti  -3.  t. 

A  ( i  5  )  =  A  ( 1  i  }  i-o(l) 

fa  (*5>  •=  3  (il)  -o(l) 

We  cannot  use  equations  !2)  and  (f..  !4)  directly  b 
in  Section  6.2,  x.d  ^  mist  be  too  close  to  each  o*h- 

In  the  case  q (z)  =  z'(iog  z)u  ,  we  typically  will  c 

c  ^5  =  log  log  1,  .  Suppose  that  we  tried  to  use 

%.\2,  and  '6.14)  by  setting  5^  (s)  =  5(c)  ,  i2^c) 

Then  we  would  find  that 

it  i  e  5  ^ 1  log2  (S^S^^  =  2(log:log  ij)  1  log2i  -  » 

It  m,  'ondit  ion  (6.6)  would  be  violated. 

Instead,  we  will  choose  a  finite  sequence 

Me)  =  5,  (-:)  <  -Me)  <  ...  <  i  (c)  =  i 

•l  —  2  —  —  n 

for  which,  for  i  =  l,2,...,n-l  , 

^  log  (1  .;■*■)-  0  as  e  -  0 

1  1  +  ±  1 

't. 12)  and  (6.14)  do  apply,  so  that  our  (joal, 
r-T-f  d  i  ately. 

-  .  -ri  ty ,  we  begin  with  an  example.  Let 


'  6 . .  1 9  , 


-1  .  1  .  .  , 

:  i  =  2  J-og  ,  log  3  , 


Define  j  .  ,  *  _  ,  5  ,  5,  as  follows 

l  z  3 


(6.20) 


(a)  51  = 

(b)  5  -v  = 


(c)  5,  = 


M  (2)  .  -1,  ,  1/4 
e  exp  (log  (-:  )  ) 

.  3.-1, 

s  log  U  ) 


(d)  >5  „  = 


Here  we  are  using  the  notation 

log2x  =  (log  x)  2  ,  log  ^2*  x 

or,  more  generally. 


=  log  log  x 


(6.21) 


.  (n+1)  .  .  (n) 

log  x  =  log  log  x 


logn+^x  =  (log  x)  (lognx) 


It  is  obvious  that  5 ^,...,6^  defined  by  (6.20) 
isfy  (6.17).  To  see  that  they  satisfy  (6.18)  we  star 
putting  £  more  explicitly  on  the  right  side  of  (6.19 
Divide  (6.19)  by  5  and  apply  log. 

legs  logo  ^  +  log  +  log  ^  5  ^  =  (l+o(l))log 
(6.22)  log'  '£  =  0(1)  +  log  '  '  5  =  (l+o  (1)  )  logs  * 

s-16  =  \  log125!-1  =  !-|  +  o(l))  log  (2 '  £~ 1 


Now  substitute  '6.22':  a.-. 


efi"1  log2(52i^L) 


=  (2-ro  (i)  )  (loij U)  -__1> [logic  i52)  ~  logic  ~-:J 
=  12«  11),  < log' 121  r1)'1  [  (log(2)c'1)L/4  - 


0 ID  -  log(3)c-1]2 


(2)  -1-1 


(2)  -Ll/4,2 


=  (240(1))  (log  c  )  [  (log'  c  ) 


=  (240(1))  (log'  t  ) 


loa(2)c-V1/2  ^0  as  e  Q 


Thus  (6.18)  holds  for  5  and  52* 

Next  substitute  (6.20  fb)  ,  (c ))  into  (6.18i 

e521log2(53521)  =  ex?(-(log(2)c_1)1/4]  [logic'1!  3) -log u"1*^)  ]2 

.  (2)  -1,1/4,  ,  (2).- 1.  n  (2)  -I  * 

=  exp  [-(log  c  )  ]  1 3  uag  (c  )-(log  c  ) 


=  (34o(1))  exp[- (log  ^  c  1)  log  ^  e  1  0 


as  c  -  0 


Thus  (6.19)  holds  for  i  ,  and  6  . 

Last,  substitute  (6.20  (c),  (d))  into  (6.18). 

eS^log2^'1)  =  log' 3  (s'1)  [log  c"1  -  3  log^’e"1]2 

=  (l+o(l))  log'1  c-1  -  0  as  c  -  0 


Thus  conditions  (6.17)  and  (6.18)  hold  in  their  entirety 
so  that  (6.15)  holds  for  this  example  of  primary  interest, 


6.4.  ac  now  give  another  example ,  one  for  which  we 
cannot  prove  (6.15) .  It  indicates  that  we  need  a  new 
assumption  in  order  to  prove  (6.15) ,  and  it  suggests  what 
form  that  assumption  should  take — a  slight  surer,  the  nine  of 
assumption  (A) . 

We  shall  construct  a  function  f  (: )  which  grows  to 
infinity  as  5  -*■  0  more  slowly  than  any  iterated  logarithm 
log  ^  5  It  is  conceivable  that  some  0($ )  =  5  *  (y  + 

+  o  (1) )  satisfying  (A)  might  have  an  error  term  so  large 
that,  to  satisfy  (3.21)  ,  we  must  have  c  ^5  <_  f  (0  )  or 

f(e).  We  shall  show  that,  for  such  a  case,  the  method  of 
subsection  6.3  will  not  work,  so  we  cannot  prove  (6.15)  for 

First,  we  construct  f  {5 ) .  For  N  =  1,2,3,...  define 


f  ^  (5  )  =  log  6 


(6.23) 


N+l 


(5)  = 


r  fN(«>  if  fN  C6)  <  N 

l  log  fN  (6 )  if  fN(5)  >  N 


For  any  given  (  (0,1),  log  fN<5)  <  fN  (3 )  ,  so  there 

is  a  finite  N  =  N(c)  for  which 


fN-l(5)  > 


W5)  = 


Therefore  we  may  define  a  limiting  function,  shown  in 
Figure  6.1. 


M 


(6.24)  f.(5)  =  lim  fN(5)  =  fN(6)(o)  =  log(N(6),5_1 

N-<-® 

Note  that  N  (& )  increases  to  infinity  as  5  -*■  0,  and  that 

by  (6.23),  f  (6)  <  N  (6 )  .  Equation  (6.23)  also  gives  a 

lower  bound  to  f  ,  for  if  N  =  N  (5 )  then  f =  log  fM  .  , 

so  that  f„  .  >  N-l.  Therefore,  fM  >  log(N-l).  In 
U~  1  N 

summary , 

(6.25)  log  (N  (5 )  -  1)  <  f_  (<S )  £  N  (i  )  . 

We  make  f  monotonic  by  defining 

(6.26)  f  (6 )  =  sup  {  f  (5')  :  5  <  5  1  <  1  } 

ao 

Because  f  (<5 )  -*•  ®  as  6  -*-0,  the  same  holds  true  for 

30 

f  (5 )  •  How  fast  does  f  (5  )  grow?  Let  5  <_  5  '  <  1.  Then 
by  (6.25), 


(6.27) 


f(o’)  <  N(6')  <  N(S)  <  1  +  ext  f  (5  )  = 

*  00 


=  1  +  log(N<5)  -  1J  r1 


Therefcre 


(6.28)  f  (6 )  <.  1  +  log(N(6)_1)  o'1 

For  any  fixed  N,  (6.2  8)  and  the  fact  that  N  (5 )  *►  a  imply 


that 


(6.29) 


as 


(s: 


<  < 


log(N> 


- 1 


5  -  0  . 


With  this  very  slowly  growing  f  (: )  ,  we  define  5  ( e )  ; 
first  define  e  (5  )  by 

(6.30)  £_16  =  f  (5 ) 

and  then  invert  to  get  o  (z)  (defined  for  all  snail  ex¬ 
cept  for  the  jump  discontinuities  at  5  ^e  =  1/f  =  1, 

■j#  1/3,  ...).  This  5(e)  is  extremely  close  to  e,  but 
it  still  satisfies  (2.14);  that  is,  z  <<  5  <<  1. 

There  is  now  so  much  distance  between  '  (e)  and  1  that 
we  can't  connect  across  that  distance,  using  the  method 
of  subsection  6.3.  Suppose  that  we  try  to  do  so,  i.e., 
that  we  construct  a  sequence  like  (6.17)  and  (6.18)  and 
try  to  use  it  to  make  a  connection.  We  shall  show  that 
our  task  would  never  end;  that  however  large  N  may  be, 
$N(e)  <  <  1  as  z  ■*  0. 

Specifically,  we  take  a  sequence  6^,  52'  «S  ^  ^  *  *  * 
which  satisfy  (6.17)  and  (6.18)  or  even  a  weakened  form 
of  (6.18),  that  is , 

(6.31)  5(e)  =  5.(e)  <  6-,(e)  <  6(e)  <  ... 

X  J 

(6.32)  e  6"1  log(6i+1  s'1)  -*■  0  as  e  -*•  0  . 


I 


Then  we  shall 


show  by  induction  that  the  following  two 
equations  hold,  for  all  i: 

(6.33)  5  .  -*■  0  as  e  -*  0 

i 

(6.34;  c_15;  <<  log^’f-1  as  e-0,  for  each  fixed  N . 

x  X 


equation  f 6 . 3  3)  is  the  result  we  want;  (6.34)  os  a  nece 
sary  intermediate  step  in  the  induction;  it  says  that  :■ 
can't  get  vary  far  from  e. 

Equations  (6.29)  and  (6.30)  taka  care  of  the  case 
i  =  1.  If  (6.33)  and  (6.34)  hold  for  seme  i,  then  by 
(6  .  32) 

(6.35)  log(ii1  5i+1>  -  olhe  1) 

Exponentiate  this  equation  to  get: 

(6.36)  oi+1  <  5i  exp  1) 

By  (6.33)  and  (6.34)  with  N  =  2  this  implies  that 

(6.37)  <5i+i  <  <51  exp  (log  (2)  5i1)  =  ^  log  6±l  -  0  . 

Now  that  we've  shown  (6.33)  for  (i+1) ,  equation  (6.34) 


remains . 


Choose  any  integer  M  >_  1,  and  require  -i  to  be 
small  that 


,  (M+3)  .-1  . 

log  5  ^  >1 


By  f  6 .  36)  and  (6.34j  with  (M+2)  and  (M  +  3)  in  place 


N, 


(6.38)  r-~ltSi+l  4  s"1*!  exP  (e ~  1<3  i >  < 


,,  (M+2)  .  -1.  ,,  (M+3 )  .  - 1 1 

<  (log  Qi  )  exp  (log  '3  x  )  = 


..  (M+2).-l  .2  .  -M+l)  -1 

(log  3^  )  <  log  3^ 


By  (6.35)  and  (6.34)  with  N  =  1  . 


(6.39)  log  3  iog  +  =  (1+0(1))  log  5 ^ 


Apply  log ^  to  equation  (6.39) 


(6.40)  log^^"^  =  (1  +  o(l)):.og(M)i~1  >  ~  log(M)  d"1 


.  (M+l)  .-1 

>>  log  5, 


Combining  (6.38)  and  (6.40)  give:;  the  result  we  had  sought 


-1  f  -  (M)  .-1 

e  6i+l  <<:  lo<?  0 


i+1 


(6.41) 


6.5.  Subsection  6.4  indicates  that  we  will  not  Li¬ 
able  to  solve  the  connection  problem  fully  if  (  j  be 

extraordinarily  close  to  .  To  give  room  for  ;  ( •: )  /  ■ 
to  grow  at  a  reasonable  rate,  we  must  have  some  control 
over  the  error  term  in  assumption  (A) .  So  we  strengthen 
(A) ,  requiring  the  error  term  to  approach  zero  at  least 
as  fast  as  some  iterated  logarithm: 

(A")  *U)  =  r1(Y  +  <>(log  (n)  i  5  |  “1)_1) 

uniformly  in  S  ,  as  |-|  -  0,  where 

n  is  a  fixed  positive  inteaer  and 

1 

“  *  <  Y  <  2  * 

This  assumption  still  covers  practically  all  familiar  ex¬ 
amples;  when  q(z)  =  zvloguz  we  may  take  n  =  1  in  (A''). 

Provided  that  (A")  holds,  we  can  define  : ( e )  by: 

(6.42)  o  —  e  log  s 

Then  for  n  >_  1 , 

(6.43)  log  5  =  (1  +  o  (1)  log  e 

It  is  easy  to  see  that  this  satisfies  (2.14).  It  also 
satisfies  all  the  other  requirements  that  we  have  put  on 
5(e)  in  earlier  sections,  for  by  (A"),  (6.43)  and  (6.42) 


t 


(6.44)  g  (6 )  =  S(log(n)  r1)-1  =  fi(log  (a)  e  1)"i 

=  o(log(n+1)rV2  =  o  (exp  (-25 /e)  )  . 


We  will  solve  the  connection  problem  by  finding 
<S  ,  5  2  /  •  •  •  satisfying  (6.17)  and  (6.18).  The  choice 
of  6 .  will  be  essentially  the  same  as  in  the  example 
in  6.3.  Set  5  =  d(e).  For  i  =  l,...,n+l,  set 


(a) 

(6.45) 

(b) 


,  (n+3-i)  -1.1/4 

5 2 =  £  exp  (log  c  ) 


t  ...  (n+2-i)  -1.3 

5  21+1  =  e  ®(1°^  £  ) 


(c) 


52n+4 


1 


We  need  only  verify  that  (6.18)  holds.  For  the  first 
step,  take  i  =  1. 

(6.46)  log2  ( 5  2  5 1  ^" )  =  e  6  1  [  log  ( c  5  1 )  log  ( 1)  1 2 

=  (log^rV1  [-log^’r1  + 

+  (log  ^^r1)  1/4  ]2 

-  (1+0(1))  (log (n+2)  -•~1)  ~1/2  -  0 


1 


By  almost  the  sate  reasoning,  we  can  verify  .'6. 
for  all  odd  5^,  except  j  =  2n+3: 


,,  .-1  ,  2  .  „  .-1  . 
C32i-rl  log  (:>  2i+2"  21+1^ 


=  (logln+2_i)c;"i)“3  [-3  log(n+3  i)E'1  + 


+  (log  (n+2-i)  e_1)  x/4]2 


=  (140(1))  (log  (n+2  x)e  ij-5/2  ^  o 


Next  we  cake  care  of  (6.18)  for  even  s j : 

(6.48)  £d  “  l  log2  (5  2i+15  2i>  =  exp[-(log(n+3":L)£"1)  i/4]  [. 


n  (n+3-i)  -1.1/4  _  .  (n+3-i)  -1,2 

...  -  (log  e  )  +3  log  t  ] 


[(3  +  o  (1)  )x]2  exp(--x^4)  -*■  0 


where  =  log^n+3  ^  e  x 


The  last  step  is  <5^  =  5 2 n + 3  »  i.o.,  i  =  n+1. 


cS2M3  ’•og2|52n+452i+3l  *  <1°^' V3 [-log  t -3  log 'V1]2 


(l+o(l))  (log  r1)-1 


0 


Therefore,  (6.18)  holds,  so  (6.15)  holds  and  the  connec¬ 


tion  is  complete. 


t 
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of  turning  points,  which  includes  logarithmic  branch  points  of  q(z),  among  many 
others.  To  this  end,  a  delicate  contraction  for  an  integral  equation  differing 
from  those  of  Langer  and  Olver  is  used  to  show  that  Bessel  functions  can  still 
1  approximate  the  solutions  at  a  certain,  small  distance  from  the  irregular  point 

of  (1.1),  even  though  not  uniformly  near  it.  A  novel  feature  of  the  analysis  is 
that  the  extreme  variation  of  the  exponential  kernel  is  here  controlled  even  on 
non-progressive  paths.  Connection  is  completed  radially  by  means  of  the  same 
integral  equation. 


